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Abstract 

We study theoretical and computational properties of the pressure function for 
subshifts of finite type on the integer lattice Z'*, multidimensional SOFT, which are 
called Potts models in mathematical physics. We show that the pressure is Lipschitz 
and convex. We use the properties of convex functions in several variables to show 
rigorously that the phase transitions of the first order correspond exactly to the points 
where the pressure is not differentiable. We give computable upper and lower bounds 
for the pressure, which can be arbitrary close to the values of the pressure given a 
sufficient computational power. We apply our numerical methods to confirm Baxter's 
heuristic computations for two dimensional monomer-dimer model, and to compute the 
pressure and the density entropy as functions of two variables for the two dimensional 
monomer-dimer model. 

2000 Mathematics Subject Classification: 05A16, 28D20, 37M25, 82B20, 82B26. 
Keywords: Pressure, density entropy, multidimensional subshifts of finite type, transfer 
matrix, first order phase transition, monomer-dimer model. 

1 Introduction 

The most celebrated models in statistical mechanics are the Ising models, introduced by 
Ising in [18], and their generalizations to Potts models [25]. Usually, the one-dimensional 
Ising or Potts models admit a closed-form analytical solution and do not exhibit the phase 
transition phenomenon, as in the case of the original work of Ising for ferromagnetism. The 
importance of Ising models was demonstrated by Onsager's closed-form solution for the 
two-dimensional ferromagnetism model in the zero-field case [23] , which does exhibit phase 
transition at exactly one temperature. Unfortunately, there are only a handful of known 
closed-form solutions for two-dimensional Potts models, including the dimer problem due to 
Fisher, Kasteleyn and Temperley [8], [20], [29]; residual entropy of square ice by Lieb [22]; 
hard hexagons by Baxter [5]. See also [6]. 

Thus, most of the interesting Potts models, in particular all problems in dimension 3 and 
up, are treated by ad hoc asymptotic expansions or by some kind of numerical solutions, 
in many circumstances with the help of Monte Carlo simulations, which usually have a 
heuristic basis. The aim of this paper is to introduce a new mathematical foundation to 
this subject, which also gives rise to reliable numerical methods, using converging upper 
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and lower bounds, for computing the pressure and its derivatives for known quantities in 
statistical mechanics. In principle, these quantities can be computed to any accuracy given 
sufficient computing power. The first-order phase transition is manifested by a jump in a 
corresponding directional derivative of the pressure, which can be detected within the given 
precision of the computation. Our approach to the phase transition is significantly simpler 
than the approaches using the Gibbs equilibrium measures corresponding to the pressure, 
e.g. [1, 7, 19, 27]. In models with one variable, the situation is relatively well understood 
by physicists. The basic argument of phase transition in Ising model is due to Peierls [24]. 
For more modern account of the physicist's approach see [14, pp' 59]. 

We now introduce the main ideas of this paper as nontechnically as possible. Assume 
that we have a standard lattice Z'', consisting of points in d-dimensional space R'' with 
integer coordinates, which we call sites. Each site i = {ii, . . . ,^£^)''" is occupied by exactly 
one particle, or color, out of the set (n) := {1, . . . , n} of n distinct colors (if we do not insist 
that every site be occupied, we agree to use color n for an unoccupied site). In general, one 
has a local type of restriction on the allowed configurations of the colors, which is called a 
subshift of finite type, or SOFT, known as the hard-core model in physics terminology. The 
exact definition of a SOFT is given in the next section. For an example of SOFT, consider 
the residual entropy of square ice studied in [22]. This entropy is the exponential growth 
rate of the number of colorings of increasing sequences of squares in with n = 3 colors, 
subject to the local restriction that no two adjacent sites receive the same color. More 
generally, we consider a nonempty near neighbor SOFT (NNSOFT), specified by a d-tuple 
r = (Fi, . . . , Td), where each F^ C (n) x (n) is a digraph whose set of vertices is the set (n) 
of colors. Two adjacent sites i and i + e^, where e/c = {6ik, ■ ■ ■ , Sdk)^ , are allowed to receive 
the colors p and q respectively only if (p, q) e Ffc. We denote the set of all allowed colorings 
in this NNSOFT by Cr(Z''). 

We assume for simplicity of the exposition that the Hamiltonian of a particle of color i 
is Ui e M. If this is not the for the Ising model or the monomer-dimer model, there 

is a way to reduce such a model to our model by enlarging the number of colors. We show 
how to carry out this reduction for the monomer-dimer model. 

For m = (mi, . . . ,mrf) e N"^, let (m) denote the rf-dimensional box (mi) x • • • x (m^). 
Let (j) : (m) — > (n) be a coloring (m) with n colors, i.e., an ensemble of vol(m) := nii ■ ■ ■ 
particles of n kinds occupying the sites in (m). Let Ci((/)) be the number of sites in (m) 
colored with color i. Let c(^) = (ci (<;/)), . . . , c„(^))''" and u = (ui, . . . , «„)''" G M". Then the 
Hamiltonian of the system (p is equal to c{(f))^u. The grand partition function corresponding 
to the set Cr((ni)) of all colorings (p ■ (m) (n) allowed by T is given by 

Zr(m,u):= ^ e'^^'^)^". (1.1) 

0eCr((m)) 

It is well-known that logZr(m, u) is a convex function. Furthermore, the multisequence 
logZr(m, u), m e N"^ is subadditive in each coordinate of m. Hence the following limit 
exists 

D / ^ r logZr(rn,u) 

Pr(u) := hm -— , (1.2) 

m^oc vol(m) 

where m — > oo means m,j — > oo for all j G (d). This limit is called the pressure function. 
The value hr ■= -Pr(O) is the (free) entropy of the corresponding SOFT, and our previous 
paper [12] was devoted to the theory of its computation. The function Pr(-) : R" — >■ R is 
a Lipschitz convex function. Hence it is continuous and subdifferentiable everywhere, and 
differentiable almost everywhere. Assume that Pr{') is differentiable at u with gradient 
vector p(u) = (pi(u), . . . ,p„(u))^. Then p(u) is a probability vector, and Pi{u) is the 
relative frequency, or proportion, of color i corresponding to the pressure Pr(u). We show 
that the points u where Pr(u) is not differentiable correspond to phase transitions of the 
first order, i.e., these are points u where the proportions of the colors are not unique. 

Let H„ C R" denote the simplex of probability vectors. Assume that p € H„ is a limiting 
color proportion vector for some multisequence of configurations in Cr((m)), m — >^ oo. Then 
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one can define the density entropy ft.r(p) as the maximal exponential growth rate of the 
number of configurations, the maximum being taken over all multisequences whose color 
proportion vector tends to p. (See for example [15] for the special case of the monomer-dimer 
configurations.) We denote by Ilr C n„ the compact set of all limiting color proportion 
vectors. 

Let Pr(-) : R" — >■ [— /ir,oo] be the conjugate of Pt{-), which is called the Legendre- 
Fenchel transform in the case of differentiable convex functions [2, 26]. Recall that P^{-) 
is a convex function. We show that for a limiting color proportion vector p, that is also a 
subgradient of the pressure function somewhere, /ir(p) = — -Pr(p)- Thus /ir is a concave 
function on each convex set of such vectors p in Ilr • 

We next show that in many SOFT arising in physical models, the set Ilr is convex and 
the function hr ■ Hr IR+ is concave. A simple example is as follows. Assume that our 
SOFT given by T has a friendly color, say n. This is, in each digraph the vertices n 
and i connected in both directions, i.e. (n, i), (i, n) € F^, for i = 1, . . . , n and k = 1, . . . ,d. 
Then Ilr is convex and /irjllr is concave. The hard core model has a friendly color. The 
monomer-dimer model has essentially a friendly color, which corresponds to the dimer, hence 
Hr = Ilrf+i and /ir|IId_|_i is concave. These results can be viewed as generalizations of the 
result of Hammersley [15]. 

For numerical computations of the pressure one needs to have lower bounds for the 
pressure, which converge to the pressure in the limit. (The convergent upper bounds are 
given by '"^yoHm)'"^ ' since the multisequence logZr(m, u), m e N** is subadditive in each 
coordinate of m.) We extend the results in [9, 12] to give lower convergent bounds if at 
least d — 1 digraphs out of Fi, . . . , F^ are symmetric. (A digraph F C (n) x (n) is called 
symmetric, (reversible), if the diedge (j,j) is in F whenever (j, «) is in F.) This condition 
holds for most of the known physical models. In this paper we show how to apply the 
computational methods developed in [12] to the pressure. We demonstrate the applications 
of our methods to tlie two dimensional monomer-dimer model on Z^. First we confirm the 
heuristic computations of Baxter [4]. Second, we find numerically a number of values of the 
pressure function P2{vi,V2) and the value of the density entropy ft.2(pi,P2) for dimers with 
densities pi , p2 in the directions xi , X2 respectively. In Figures 1 and 2 we give the plots of 
these functions. These computations go beyond the known computations of [4, 16], where 
one considers the total density of dimers, (which reduce to the computations of functions of 
one variable). 

We hope to show that the methods of this paper can be applied to other interesting 
models. We already know that our approach works for the numerical computation of the 
pressure function for the 2D and 3D Ising models in external magnetic field. (It is similar 
to the monomer-dimer model we study here.) We plan to study if our numerical compu- 
tations are precise enough to discover the second order phase transitions, which occurs in 
multidimensional Ising models. 

We now survey briefly the contents of the paper. In Section 2 we describe in details 
SOFT, NNSOFT and the pressure function. We also show that in the one-dimensional case, 
the pressure is the logarithm of the spectral radius of a corresponding nonnegative matrix. In 
Section 3 we show that under certain symmetry (reversibility) assumptions on d— 1 digraphs 
among Fi, . . . , F^, we have computable converging upper and lower bounds for the pressure. 
In Section 4 we relate certain properties of a convex function, as differentiability and its 
conjugate Prip), (the Legendre-Fenchel transform), to the physical quantities associated 
with a given SOFT, i.e. the corresponding Potts model. In particular we show that the 
points where the pressure Pr{-) is differentiable correspond to unique color frequency vectors. 
On the other hand the points were the pressure is not differentiable correspond to the phase 
transition of flrst order, since to this value of u correspond at least two different color 
frequencies. We also relate the density entropy /ir(p) to the conjugate function Pr(p)- 
Section 5 we apply the results of Section 4 to a one-dimensional SOFT. The importance 
of one-dimensional SOFT is due to the fact that our approximations of the pressure are 
obtained by using the exact results on one-dimensional SOFT. Section 6 we apply some of 
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our results in Section 4 to the monomer-dinier model in H^. We also relate our results to the 
works of Hammersley and Baxter [15, 4]. This is done by using the fact that the monomer- 
dimer model in can be realized as SOFT with 2d + 1 colors [10, 12]. As we pointed 
out in [12] this SOFT docs not have symmetric properties, and hence can not be used for 
computation. In Section 7 we use the symmetric encoding of the monomer-dimer model 
developed in [12], to obtain computer upper and lower bounds for the pressure function. 
In Section 8 we apply our techniques to the computations of two dimensional pressure and 
density entropy for the monomer-dimer model in 1? . 

2 SOFT, NNSOFT and Pressure 

We use the notation (r) := {1, . . . , r} for r € N := {1, 2, 3, . . .}, and for m = (mi, . . . , md) S 
N'', (m) := (mi) x • • • x {md) denotes a box with volume vol(m) := mi • • • rud- Then 
(n) is the set of all colorings (f> : (m) (n) of (m) with colors from (n) . We denote by 
c{(j))i := the number of sites in (m) colored with the color i S (n), and let c(^) := 

{c{(j))i, . . . ,c{(f>)n)^ . Similarly, with Z := {0,±1,±2, . . .}, (n)^ is the set of all colorings 
: Z'^ — >• (n) of Z'^ with colors from (n). Given a d-digraph T = (Fi, . . . , F^) on (n) x (n), 
let Ct{'^'^) C (n)^ be the set of all T -colorings^ namely colorings (j) = (</'m)mez<' '= i^)^ 
such that for each i & Z'^ and fee (d), ((/>i, (/>i+efc) S F^, where e^, is the unit vector with 
fcth component equal to 1. In ergodic theory, the set Cr(Z'') is called a nearest neighbor 
subshift of finite type (NNSOFT). 

A general SOFT can be described as follows. Let M e N'' and a nonempty subset 
V C (n) be given. Every element a gP is viewed as an allowed coloring (configuration) 
of the box (M) with n colors. For i € Z'^, we define the shifted coloring ri(a) of a (E "P as 
the coloring of the shifted box (M) + i that gives to the site x + i the same color that a 
gives to X e (M). We denote by Ti{V) the set {Ti(a) : a G V}, and regard it as the set of 
allowed colorings of (M) + i. A coloring cj) G (n)^ is called a V-state if for each i G Z'' the 
restriction of (j) to (M) + i is in Ti{V). We denote by {n)^'' (V) the set of all P-states. In 
ergodic theory, the set (n)^ (P) is called a subshift of finite type (SOFT) [28]. 

Each NNSOFT Cr(Z'') is a special kind of SOFT obtained by letting M (2, . . . , 2) 
and V the set of all colorings (j) G (n)^"^^ such that i. i + G (M) imply (^i, (jji+ek) € F^. 
Conversely [9], each SOFT {n)^\v) can be encoded as an NNSOFT C'riZ'^), where T = 
(Fi, . . . , Frf) is defined as follows. Take N = and use a bijection between V and {N). 
The digraph F^ C [N) x {N) is defined so that for a,b G P we have (a, b) G F^ if and only 
if there is a configuration cj) G (n)^'^+''''^ such that the restriction of (j) to (M) is a and the 
restriction of (f) to (M) + Cfc is Te^{b). Because of this equivalence, we will be dealing here 
with NNSOFT only. 

In the sequel we will take limsup and liminf of real multisequences (am)meN'' as m — > oo. 
In order to be clear, we define these here and observe that they are limits of subsequences 
[12]. We also define the limit of real multisequence in terms of limsup and liminf, which is 
equivalent to other definitions in the literature. 

Definition 2.1 Let (am)meN'' multisequence of real numbers. Then 

(a) lim supjj,_j.oQ Om is defined as the supremum (possibly ±qo) of all numbers of the form 
limsupg_^oo flm, > where (mg)qgN is a sequence in satisfying Viuvq^^ = oo, i.e., 
limq^oo(mg)fe = oo for each k G (d). We define lim inf ^^-cxa oim similarly. 

(b) limm^^oo (im = a means limsup^j.^^ am = liminf m^-oo am = a. 

As in [12], given an NNSOFT Cr(Z'0 and m G N'', we denote by Cr((m)) the set of all 
colorings </> G (n)^™^ such that 1,1 + 6^ G (m) imply ((/>!, (/)i+efc) G F^. Similarly, we denote 
by Cr,top((m)) C Cr((m)) the projection of Cr(Z'') on (m), i.e., the set of colorings in 
(n)<""> that can be extended to colorings in Cr{Z'^), and by Cr,per((m)) C Cr,top((m)) the 
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Zr,per(m,u) := ^ 



set of periodic F-colorings with period m, i.e., the set of colorings in (n)^™^ that can be 
extended to colorings in Cr(Z'') with period m. For a weight vector u = {ui, . . . , G M" 
on the colors, we define 

Zr(m,u):= Yl e'^*'""' (^.l) 

<^eCr({m» 

Zr,top(m,u):= e'^^^)"", 

06Cr,top(<m>) 

pc(0)^u 
,^eCr,per((m)) 

As usual, a summation over an empty set is understood as 0. Obviously 

#Cr((m)) = Zr(m,0), 
#Cr,top((m)) = Zr,top(m,0), 
(m, 0). 

A function /(u) > on R" is called log-convex when log /(u) is convex. (The zero function 
is by definition log-convex.) Recall that the log-convex functions are closed under linear 
combinations with nonnegative coefficients [21]. Since for c G R" the function e*^ " is log- 
convex, the weighted sums Zr(m, u), Zr,top(m, u), Zr,por(m, u) are log-convex functions of 
u for each m e N''. As in [12] it follows that for a fixed u, log.Zr(m, u) and log.^r,top(m) u) 
are subadditive in each coordinate of m, and so the limits (2.2) and (2.3) below exist. The 
quantities 

Pr(u) := lim l^i^L^L^, (2.2) 
nn-oo vol(m) 

log2r,top(m,u) 



Pr,top(u):= lim ^' ^V ' ^^ (2.3) 

itn-oo vol(m) 

D i \ V logZr,per(m,u) 

-Pr,per(u) := limsup '-^--^ (2.4) 

m-)-cx> VOl(m) 

are called the pressure, the topological pressure and the periodic pressure of Cr(Z''), re- 
spectively. The special cases hr := -Pr(O), /ir.top := -Pr,top(0), /ir.per := -Pr,per(0) are the 
entropy, the topological entropy and the periodic entropy, respectively, discussed in [12]. 
Clearly 

-oo < Pr,pcr(u) < Pr,top(u) < Pr(u). 

By the log-convexity of Zr{m, u), .Zr,top(ni, u), Zr,per(m, u), it follows that -Pr(u), Pr,top(u), 
-Pr.pcr(u) arc convex functions on K". (We agree that the constant function — oo is convex.) 
As in [9], one has the equality 

Pr,top(u) =Pr(u). 
Since logZr(m, u) is subadditive in each coordinate of m, it follows that 

< ^o,Zr im n) _ 
vol(m) 

In the one-dimensional case d = 1, we can express Pp^ (u) as the logarithm of the spectral 
radius (largest modulus of an eigenvalue) of a certain n x n matrix as follows. 

Proposition 2.2 Let = idij)i,je{n) be the adjacency matrix ofVi, and let 

^ri(u) = (rfjj(u))ije(„) he defined by 

dij{u) :=di, -e^(^^"+^7"). (2.6) 
Let pri(-D(u)) be the spectral radius of Dri{u). Then 

Pr,(u)=logpri(£>(u)). 
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Proof Recall the following characterization of the spectral radius of a nonnegative matrix 
M: for any vector w with positive components, we have p(M) = limfe_>.(x,(w^M'°w)fc (cf., 
Proposition 10.1 of [10]). Consider the positive vector 

l(u) = (eK""),e(„). 

Since Cri((mi)) is the set of walks of length mi — 1 on Fi, we have Zr^{{mi),u) = 
l(u)"^£>ri(u)'"i-il(u). Therefore 



logZri(mi,u) 
hm ^ = Pn (u) ■ 

n^-oo mi 



logpri(-D(u)) = lim 

mi— ^OO 



= lim 

mi—>oo 



logl(u)'Dr,(uril(u) 



mi 



logl(u)'^£)ri(u)™i-il(u) 



mi 



3 Main Inequalities for Symmetric NNSOFT 

In this section we derive bounds for the pressure analogous to those for the entropy in [12, 
Section 3] under the assumption that some of the digraphs Pi, . . . , Fd arc symmetric. 

For d > 2, consider m = (mi,...,md) S N'' and m~ := (m2,...,md)- We denote 
by T(m) the discrete torus with sides of length mi, . . . ,md, i.e., direct product of cycles 
of lengths mi, . . . ,ind- Let Cr,pcr,{i}(m) be the set of F-colorings of the box (m) that 
correspond to F-colorings of T(mi) x (m~), i.e., that can be extended periodically in the 
direction of ei with period mi into F-colorings of Z x (m~). We can view these colorings 
as F-colorings of the box (m~), where F = (F2, . . . , P^), for each k the vertex set of P^. is 
the set P™per of closed walks a = (ai, . . . ,ami,ai) of length mi on Pi, and (a, 6) e Pfe if 
and only if (oj, bi) e Pfe for i = 1, . . . ,mi. For this reason, the limit (3.1) below exists and 
is equal to the pressure (u) of the NNSOFT Cf{Z'^-'^): 

Zr,per,{i}(m,u):= ^ e-^f"^)"", 

<)!'eC'r,per,{i}(m) 

D / ^ r log^r,per,{i}(m,u) 

Pr(mi,u):= lim \ , mi e N. (3.1) 

m-^-oo vol(m~) 

Then Pr(mi,u) is a convex function of u G R". In the degenerate case mi = 0, we 

define Zr,per,{i}((0, m~), u) to be #Cp(m~) (regardless of u), where Cp(m~) is the set 

of (F2, . . . ,Fd)-colorings of the box (m~). Then (3.1) is also valid for mi = 0, where 
Pr(0,u) := P(r„...,r.)(0) is the entropy of C(r„...,r,)(Z''-'). 

Theorem 3.1 Consider the NNSOFT Cr(Z^) for d>2, and let Pr(u) and Pr(rni,u) 
be defined by (2.2) and (3.1), respectively. Assume that Fi is symmetric. Then for all 
p, r e N and g S Z+, 

J'r(2r,u) Pr(p + 2g,u)-Pr(2g,u) 

> Pr(u) > . (3.2) 

Proof Fix m~ = (m2, . . . ,md) e N"^"^ and let f2i(m~) be the following transfer digraph 
on the vertex set Cp(m"), analogous to the transfer digraph Vld^m.') described in [12, 
Section 1]. Vertices i and j satisfy (i, j) e Oi(m~) if and only if [i, j] G Cr(2, m~), where [i, j] 
is the configuration consisting of i, j occupying the levels Xi = 1,2 of ((2, m~)), respectively. 
Let N = #(7p (m~) and let Z)r(m~) = (^ij)i j^cf (m-) be the NxN {0, l)-adjacency matrix 
of f2i(m-). Let £)r(m-,u) = (rfij(u))i j£c?-(m-) be defined by 

rfij(u) = dii ■ e"*" i,j e Cp (m-), (3.3) 
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and let the positive vector l(u) be defined by 

l(u) = (e^)igc-(m-)- 

Then 

l(u)TDr(m-,u)™il(u) = Zr((mi,m-),u), 
and as in the proof of Proposition 2.2 

iogp{Dr(m ,u))= lim — ^ 

mi->C!0 mi 

(In particular, p{Dr{m^ , u)) is a log-convex function of u [21].) It follows that 

logp(Dr(m-,u)) ^ ^.^ logZr((mi, m~), u) 
vol(m~) mi->oo mivol(m~) 

Now send TO2, . . . , ma to oo, and observe that by (2.2) and (2.5), the right-hand side of (3.4) 
converges to fr(u) and bounds it from above for each m~. Thus we obtain an analog of [9] 

logp(IJr(m-,u)) ^ ^_ ^ 

vol(m~) 
lim logp(Z)r(m-,u)) ^^ 

m-^-oo VOl(m~) 

Next, we observe that 

tr£)r(m-,u)« = Zr,per,{i}((«,m-),u), q€l+, (3.7) 
where £)r(m~,u)° is the N x N identity matrix. Recall that 

N 



tr£)r(m-,u)'' = 5^A«, qeZ+, 



i=l 

where Ai,...,Ajv be the eigenvalues of Dr(m~,u). Since £>r(m~,u) is a nonnegative 
matrix, its spectral radius p{Dr{m~ ,uj) := max^g^^^ |Ai| is one of the A^ by the Perron- 
Probenius theorem. Since by assumption Fi is symmetric, rii(m^) and hence £'r(m^.u) 
are symmetric. Therefore Ai, . . . , Ajv are real, and hence tr£'r(in~, u)^'' > p(I>r(m~, u))^'' 
for each r e N. Taking logarithms and using (3.7), we obtain 

log^r,per,{i}((2r, m-),u) ^ logp(Jr(m^, u)) 
2rvol(m~) ~ vol(m~) ' 

Sending TO2, • • • j md to oo in (3.8) and using (3.1) and (3.6), we deduce the upper bound for 
Pr(u) in (3.2). 

To prove the lower bound in (3.2), we note that 

trDr(m-,u)^+^« = E ^r'" < El^^l'^'^ = El^^l'^?' 

i i i 

< Ep(^r(m-,u))fA^« =p(£>r(m-,u))f trDr(m-,u)29 

i 

and thus by (3.7) 

f -A tri^r(m-,u)P+2-^ ^r,pcr,{i} ((p + 2g, m'), u) 

^ tri;,(m-,u)^. - Z,,p.,,,,((2,,m-),u) " ^'''^ 
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Therefore 



logp(Dr(m ,u)) ^ log^r,per,{i}((p + 2g,m ), u) - log Zr,per,{i}((2g, m ),u) 
vol(m~) ~ pvol(m~) 

Sending to oo and using (3.6) and (3.1) (recall that the latter holds for mi G Z+), we 
deduce the lower bound in (3.2). ■ 

When d = 2, Pr(mi,u) is the pressure of the NNSOFT Cp^(Z) (recall that Pr(0,u) 
is the entropy hr2)- Since this is a 1-dimensional NNSOFT, Proposition 2.2 implies that 
Pr{mi,u) = log/o(Z)p^(u)), where £)p^(u) is defined as in (2.6). We denote p(£'p^(u)) by 
^2(^1, u), and obtain the following corollary to Theorem 3.1. 

Corollary 3.2 Let d = 2 and assume that Fi is symmetric. Then for all p,r gN and 

log^2(2r.u) ^ ^ ^^^^ ^ Ioge2(p + 2q,u)-log02{2q,u) 

^ > Pr(u) > , (3.1U) 

where 62 is defined above. ■ 

In (3.10) take g = and p = 2r, and send r to 00. Clearly the upper and lower bounds 
then converge to -Pr(u). Hence -Pr(u) is computable, as shown in [9] for the entropy -Pr(O). 
Combining the arguments of the proof of Theorem 3.1 with the arguments of the proof of 
Theorem 3.4 in [12], we obtain 

^ logfl(_Dr(rn^. u)) 

Pr(u) < -. r- , TO2, . . . , TOrf even, F2, . . . , F^ symmetric, (3.11) 

vol(m-) 

where £'r(m~,u) is defined in (3.3). 



4 The Conjugate of Pressure and the Density Entropy 

The purpose of this section to exhibit a striking connection between the conjugate function 
Pp(p) of the pressure Pr(u) and the density entropy /ir(p) for certain probability vectors 
P- 

First we need to recall some properties of convex functions, which can be found in [26]. 

We adopt the notations of that book. 

In this paper we consider only convex functions / : M™ — > M U {+00} that are not 
identically equal to +00. Such convex functions are called proper in [26]. Let / be a proper 
convex function. Then dom/ := {x € : /(x) < 00}, which is called the effective 
domain of /, is a nonempty convex set in M™. Let L C M™ be the minimal affinc subspace 
that contains dom /, and let I be its dimension. The affinc transformation A that maps 
L onto maps dom/ onto a convex set C C M'. We denote the interior of C by int C. 
Then ri (dom/) := yl~^(int C) is called the relative interior of dom/ (note that if / = 0, 
then dom/ and ri (dom/) consist of the same single point). A proper convex function 
/ is Lipschitzian relative to any closed bounded subset of ri (dom/) [26, Thm 10.4]. In 
particular, / is continuous on ri (dom /) [26, Thm 10.1]. 

A proper convex function / is called closed if / is lower semi-continuous [26, Section 7, p. 
52] . In particular, if / : M™ ^ M is convex, then dom / = M™ = ri (dom /) , / is a continuous 
function on M™, hence closed, and / is Lipschitzian relative to any closed bounded subset 
of R™. 

A vector y G R™ is called a subgradicnt of f at x ^ M™ if /(z) > /(x) + y^(z — x) for 
all z G M™. The set of all subgradients y at x is called the suhdifferential of f at x and is 
denoted by 9/(x). As usual, for any set 5* C R™, df{S) denotes [Jxesdf{x.). Obviously 
9/(x) is a closed convex set. If 9/(x) 7^ 0, then / is said to be subdifferentiable at x. A 
proper convex function / is not subdifferentiable at any x ^ dom/, but is subdifferentiable 
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at each x € ri (dom/) [26, Thm 23.4]. Recall that / is dijferentiable at x if there exists a 
vector V/(x) = y € (necessarily unique) such that /(x + w) = /(x) + y^w + o(||w||), 
w — >■ 0. The vector V/(x) is called the gradient of f at x. We denote by diff / the set of all 
points where / is diffcrcntiablc, so 9/(diff /) denotes the set of all gradient vectors of /. A 
proper convex function / is differentiable at a point x G dom / if and only if 9/(x) consists 
of a single point, which is then V/(x) [26, Thm 25.1]. 

Assume that / is a proper convex function and int (dom /) ^ 0. Then diff / is a dense 
subset of int (dom/), / is differentiable a.e. (almost everywhere) in int (dom/), and V/ is 
continuous on diff / [26, Thm 25.5]. Moreover, for each x e int (dom/) \ diff /, the convex 
set 9/(x), which consists of more than one point, can be reconstructed as follows from the 
values of the gradient function V/ on diff /. Let ^(x) consist of all the limits of sequences 
V/(xj) such that Xj e diff / and Xj x. Then 5(x) is a closed and bounded subset of 
and 9/(x) = conv 5(x) [26, Thm 25.6, 7.4]. 

We now recall properties of the conjugate of a convex function / [26, Section 12], denoted 
by /*: 

/*(y) := sup x^y - /(x) for each y e R™. 

Since we assumed that / is proper, it follows that /* is a proper closed convex function; 
moreover, if / is closed then /** = / [26, Thm. 12.2]. A straightforward argument shows 

that 

/*(y) = x''"y — /(x) for each subgradient y e 9/(x). (4.1) 

Recall that if / is closed then df* is the inverse of df in the sense of multivalued 
mappings, i.e., x G df*{y) if and only if y € df{x) [26, Cor 23.5.1]. In what follows we 
need the following result: 

Lemma 1 Let f be a proper closed convex function on M™. Then df{R"^) is exactly the 

set of points in R™ where f* is subdifferentiable. In particular, ri (dom/*) C df{W^) C 
dom/*, and the closure o/9/(R™) is equal to the closure of dom f*. □ 

Proof Assume that /* is subdifferentiable at y. Then there exists some x e M™ such that 
X e df*{y). Since / is closed, it now follows that y e 9/(x), so y e df{W"). The converse 
is shown in the same way. The first statement of the "in particular" now follows since /* is 
proper and hence, as noted above, is subdifferentiable in ri (dom /* ) but not outside dom /* . 
The second statement of the "in particular" follows from the first one and from the fact 
that for a convex set S such as dom/*, ri S and S have the same closure [26, Thm 6.3]. ■ 

We return to a general NNSOFT Cr(Z''). We assume throughout that Cr(Z'^) 7^ 0, for 

otherwise the pressure function Pr is identically —00. 

Proposition 4.1 The pressure Pr is a convex nonexpansive Lipschitz function on 
with constant at most 1 with respect to the norm \\{v-i, . . . , Vn)\\oo '■= ^8iXie{n} l^il •' 

|Pr(u + v)-Pr(u)| < llvlloo, u,veK". (4.2) 

In particular, Pr is finite throughout R" (this also follows from (2.5)). Therefore it is a 
proper closed convex function. 

Proof The convexity of Pr was pointed out in Section 1. Let (j) € Cr((m)). Then 

|c(0)^vl < vol(m)||v||oo. 
Therefore a term-by-term comparison gives 

e-™'(™)"^"~Zr(m,u) < Zr(m,u-hv) < e™^("")ll^ll~Zr(m, u). 
Take logarithms and divide by vol(m) to obtain 

logZr(m, u -|- v) logZr(m, u) 



vol(m) vol(m) 
Letting m — )• 00, we deduce (4.2). 



< IMIcx,. (4.3) 
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Since Pr : M" ^ M is a convex function, it is differentiable almost everywhere. Consider 
the set of probabihty distributions on the set of colors 

n„ := {p= {pi,...,Pn) > 0, pi H \-pn = 1}. 

For m e N, we denote 

n„(m) := {c = {ci, . . . ,Cn) e : ci -\ h c„ = m} = mn„ n Z". 

Let 

Cr((m), c) := {0 e Cr((m)) : c(?!>) = c}, for all c e n„(vol(m)). 
This is the set of r-colorings of (m) with color frequency vector c. 

Definition 4.2 A probability distribution p e n„ is called a density point o/Cr(Z'^) 
when there exist sequences of boxes (nig) C N** and color frequency vectors Cq G n„(vol(mq)) 
such that 

m„-)-oo, Cr((mg),Cg) ^ e N, and lim ^ = p. (4.4) 

9-^oo VOl(mq) 

We denote by Ur the set of all density points o/Cr(Z''). For p € Ilr we let 
. / X 1- log#Cr((mq),Cg) 

'ir(p) := sup limsup ^ — > 0, (4.5) 

m,,c, g^-oo VOl(mg) 

where the supremum is taken over all sequences satisfying (4-4)- One can think o//ir(p) o-s 
the entropy for color density p, called here the density entropy. 

It is straightforward to show (using a variant of the Cantor diagonal argument) that Ilr is 
a closed set. Furthermore, hr is upper semi-continuous on Ilr, because it is defined as a 
supremum. 

Theorem 4.3 Let Pp be the conjugate convex function of the pressure function Pr- 

Then 

(a) hr{p) < -Pr{p) for all p e Hr. 

(b) 

Pr(u) = max(p^u + /ir(p)) for all u e M". (4.6) 
penr 

For u e M", we denote 

nr(u) := arg max {p^u + hr{p)) = {p e Hr : Pr (u) = p"^u + /ir(p)}, (4.7) 

pellr 

(c) For each p e nr(u), hr{p) = -P^{p). 

(d) nr(u) C aPr(u), In particular, if u G diffPr, then Ur{u) = {VPr(u)}. Therefore 
(9Pr(diff Pr) C Hr. 

(e) Let u e M" \ diff Pr, and let S{u) consist of all the limits of sequences VPr(ui) such 
that Ui e diff Pr and Ui u. Then S{u) C nr(u). 

(f) conv nr(u) = conv S{u) = dPr{u). Hence 9Pr(M") C conv Hr C n„. 

(g) conv Hr = domPp. 
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Proof First we show that Pr(u) > p^u + /ir(p) for all p € Hr- Fix p € Hp and let m.g, Cg, 

g e N, be sequences satisfying (4.4). Wc have Zr(mq,u) > #Cr(mq, Cq)e'^'J since the 
right-hand side is just one term of the sum represented by left-hand side. Take logarithms, 
divide by vol(mg), take limsupg_^(^ and use the definition of fr(u) and the limit in (4.4) 

to deduce fr(u) > p^u -I- limsup^.^^^ ^"^"to^m") '""^ • ^'^^ take the supremum over all 
sequences mg,Cg satisfying (4.4) and use (4.5) to obtain 

Pr{u) > p^u + hr{p) for all p e Hp (4.8) 

and thus 

Pr(u) > sup p"^u + /ir(p). (4.9) 
penr 

On the other hand, (4.8) can be written as p^u— Pr(u) < — /ir(p) for all p e Hr, and then 
taking the supremum over u gives P}{p) < — /ir(p) < oo for all p G Hr- Thus we have 
established (a) as well as 

Hr C domP^. (4.10) 

We now show that for each u e K" there exists p(u) e Hr satisfying Pt(u) < p(u)''"u + 
/ir(p(u)), which together with (4.9) will establish (4.6). Observe first that 

#n„ (m) = f ™ + " ;^)^ Oim''- 1 ) , m ^ oo. 



n-1 



Therefore for each m s N'', 



Zr(m,u) =0(vol(m)"-i) max #Cr((m), 0)6*^^". 

cen„(voi(m)) 

Let 

C(m,u):=arg max #Cr((m), 0)6*=^". (4.11) 
cen„(voi(m)) 

Then for c(m, u) e C(m, u) we have 

Zr(m, u) = 0(vol(m)"-i)#Cr((m), c(m, u))e'=('"'")"'". (4.12) 

Since Cr{Z'^) ^ 0, for each m e N'^ and u G R", is a well-defined point in n„. 

Choose a sequence — )• 00 such that converges to some p(u). We have p(u) e Ilr 

by Definition 4.2. Apply (4.12) to rtig, and use the definition of Pr(u) and /ir(p(u)) to 
deduce 

D ^ ^ / / ^T , ,. fog#C'r((mg),c(mg,u)) . , , , t f \\ 

-Pr(u) < p(u) u + hmsup — — < p(u) u + /ir(p(u)), 

q_j.oo vol(mq) 

which is the desired inequality establishing (b). 

By the definition of Ilr (u), for each p e nr(u) we have — /ir(p) = P^u— Pr(u) < Pp(p). 
Combining this with (a), we deduce (c). 

Let p G nr(u) and v e M". Then the maximal characterization (4.6) of i-t(u + v) and 
(4.7) give 

Pr(u + v) > p'^(u + v) + /ir(p) = p"^v + Pr(u). 

This proves (d). 

We now prove (e). Assume that u <E M" \ diff Pp and p e S{u). Then there exists 
a sequence G diff Pp such that Uj ^ u and VPr(ui) p. We have {VPr(ui)} = 
nr(ui) C Hr by (d), and since Hr is closed, p G Hr- By definition of nr(ui) we have 
Pr(ui) = VPr(ui)^Uj + /ir(VPr(uO). When « ^ 00 we have firstly Pr(ui) Pr(u) by the 
continuity of Pr, secondly VPr(ui)^Ui — )■ p^u, and thirdly limsup /ir(VPr(ui) < /ir(p) 
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by the upper semi-continuity of hr- Therefore Pr(u) < p^u + hr{p)- This, the fact that 
p e Hr, and (4.6) show that fr(u) = p^u+ hr{p), which by definition means p € nr(u). 

We show the first first identity of (f). Let u e M" \ diff Pp, and let S{u) be as in (e). 
By (e) we have S{u.) C nr(u), and therefore 9Pr(u) = conv S'(u) C convnr(u). Since 
dPr{u) is convex, from the first claim of (d) we obtain dPr{u) I) conv nr(u). Hence 
dPr{u) = convnr(u). Clearly, convnr(u) C convHr- Hence aPr(K") C conv Hr- The 
second inclusion of the second claim of (f) follows from Hp C n„, which holds by definition 
of Hr. 

We finally show (g). By (4.10) and the convexity of domPp, we have conv Hp C domPp. 
It is left to show that domPp C conv Hp. By Lemma 1 9Pr(M") is the set of all points 
where Pp is subdifferentiable. In particular ri (domPp) C 3Pr(K"), so by (f) we have 
ri (domPp) C conv Hp. Apply the closure operator to both sides of this inclusion. On the 
left we get cl (domPp) by the convexity of domPp (by [26, Thm 6.3], every convex set C 
satisfies cl (ri C) = cl C). On the right we get conv Hp because Hp is closed. So we obtain 
domPp C cl (domPp) C conv Hr, as required. ■ 

As hr = Pr(0), we obtain from (4.6) the following generalization of [12, (4.12)], which 
deals with the case of monomer-dimer entropy: 

Corollary 4.4 

hr = max hr(p)- 

penr ■ 

For each p G Hr(M") := UueR"nr(u) we have hr(p) = -Pr(p) by (c) of Theorem 
(4.3). Since Pp is a convex function, we obtain the following generalization of the result 
Hammersley [15]. 

Corollary 4.5 The function hr{-) ■ Hp M+ is concave on every convex subset of 
Hr(M"). ■ 

To obtain the exact generalization of the result of Hammersley that Hp is convex and 

hr{-) ■ Hr — > is a concave function on the entire Hr, we need additional assumptions 
on the digraph F, which do hold for the T that codes the monomer-dimer tilings of Z''. For 
m e N'', if a : (m) — )• (n) is a coloring of a box (m) and j G U^, then to color the shifted 
box (m) + j by a means to give to x -|- j the color q;(x) for each x G (m). Recall that 
Cr((ni)) denotes the set of all F- allowed colorings a : (m) — )• (n), that is to say, such that 
if x,x + ei e (m), then (a(x),a(y)) e Fj. 

Definition 4.6 For a given digraph F = (Fi, . . . , F^) on the vertex set {n), a set T = 
UmGNt'C'r ((m)), where Cr((ni)) C Cr((m)) for each m G N'', is called friendly if the 
following condition holds: whenever a shifted box is cut in two and each part is colored by a 
coloring in F, then the combined coloring also belongs to T . More precisely, let m, n e N"^ 
anrf j G be such that (m) n ((n) -|-j) = 0, and such that T := mU ((n) -|-j) is a box (k) -I- i 
for some k G N'^ and i G Z''. Let a G C7r((m)), /S G C'r((n)), and let j : T (n) color (m) 
by a and (n) -|- j by (3. Then the coloring J : (k) — >• (n) defined by (5(x) = 7(x -|- i) belongs 
to Cr((k)). 

The digraph F is called friendly if there exist a friendly set T = UjngNdCr ((ni)) and a 
constant vector b G N'' such that if any box (m) is padded with an envelope of width bi in 
the direction of et, then each V -allowed coloring of (m) can be extended in the padded part 
to a coloring in T . More precisely, for each m G N'' and each a G Cr((m)), there exists a 
coloring in C'r((m -|- 2b)) that colors (m) -|- b by a. 

Example 4.7 Let F = (Fi, . . . ,Vd) be a coloring digraph with vertex set (n). Then F 
is a friendly digraph with b = 1 if one of the following conditions holds: 
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(a) r has a friendly color f € (n) , i.e., for eaeh i G (d) we have {f,j), ij,f) G for all 
j e (n) (we can take Cr(m) to be those T-allowed colorings 0/ (m) whose boundary 
points are colored with f). This example is useful for the hard-core model with n = 2 
and ri = {(l,l),(l,2),(2,l)}, / = 1. 

(b) T is the digraph associated with the monomer-dimer covering as defined in (6.1) (we 
can take Cr((ni)) to be the set of tilings of (m) by monomers and dimers, i.e., the 
coverings in which no dimer protrudes out of (m), as in Hammer sley ) . 

The following theorem strengthens Theorem 4.3 and generalizes the results of Hammer- 
sley in case F is a friendly digraph. 

Theorem 4.8 Let V = (Fi, . . . , r^) be a friendly coloring digraph. Then 

(a) Ilr is convex. Hence Hr = dom . 

(b) hr{-) ■ Hr — ^ 1R+ is concave. 

(c) For each u e R", nr(u) = aPr(u). 

(d) For each u € M", hT{-) is an affine function on dPr{u). 

(e) /ir(p) = --Pr(p) /^'^ 6*^^/1 p e Hr- 

Proof We first first prove (a). Let a e Cr((m)), let c(a) = (ci, . . . ,c„) S n„(vol(m)) be 
the color frequency vector of a, and let p := :;^o^^c{a). We assert that p e Hr. For k = 

(fci, . . . , kd) € N'^, we define k • m := (fcimi, . . . , kdmd) and view (k • m) as a box composed 
of vol(k) boxes isomorphic to (m), i.e., as (m) duplicated by a factor of k. Wc color each of 
these boxes by a, obtaining a coloring a(k • m) of (k • m). Clearly, p = voi(k.m) '^('^(^ ' ™))- 
Since a e Cr((m)), it follows that a(k • m) belongs to Cr(k • m), so in particular is F- 
allowed. Choosing a sequence k^ — > 00, wc deduce that p = limg^oo voi(k^-m) ^i'^(^Q ' ™))- 
Hence p € Ilr according to (4.4) , as asserted. 

Let 13 e (7r((n)). By the above argument we also have q := e Hr- We assert 

that all i,jGN satisfy j^P + j+jl ^ Ilr. Let a{n ■ m) and ,5(m • n) be defined as above. 
Notice that (n • m) is isomorphic to (m • n). By the above argument, a(n • m) e Cr(n • m) 
and /3(m • n) s Cr(m • n). We define k := (mini, . . . ,md-ind-i, {i + j)mdnd) and view 
the box (k) as composed of i + j boxes isomorphic to (m • n) aligned side-by-side along the 
direction of e^. Color the first i of these boxes by a(m ■ n) and the last j by /3(n • m), 
obtaining a coloring 7 of (k), which satisfies ^i^'^^''') ^ ^ i+j*^' T ^ C'r((k))- 

so in particular 7 is F-allowed. By the above argument we obtain that j^ijP + T+jQ ^ Hr, 
as asserted. Since Hr is closed we deduce that ap -|- (1 — a)q e Hr for all a e [0, 1]. 

Let Ilr be the convex hull of all points of the form ^^^^^ c(a) for some m and some a G 

Cr((m)). By the argument above we have Ilr Q Hr- Let p G Hr- By Definition 4.2 there 
exist sequences m.q 00 and color frequency vectors Cg G n„(vol(mq)) satisfying (4.4). Let 
ag e Cr{{mg),Cg). By Definition 4.6, a, can be extended to a coloring 5^ e C'r((mg -|- 
2b), Cg) for some Cg. Since b is constant and rtig 00, we have limg^oo voi(m^-i-2b) ^g ~ P" 
^^^^^ voi(m^-F2b) ^g ^ we have p e cl Hr- Thus Hr C Hr C cl Hr- Applying the 
closure operator, we deduce Hr = cl Hr, and since Hr is convex, so is Hr. The equality 
Hr = domPp follows from part (g) of Theorem 4.3. 

We now prove (b). Choose any £ > 0. Let p € Hr. By Definition 4.2 there exist 
sequences — >■ 00 and color frequency vectors £ n„(vol(mg)) satisfying (4.4). By (4.5) 
we may assume by selecting appropriate subsequences that the following limit exists and 
satisfies limg^oo '"^^TOU^m"^^ ''^^ — '^r(p) — s. Each coloring ag G Cr((mg),Cg) can be 
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extended to sonic coloring a.q G Cr((mq + 2b)). Denote by c(aq) the color frequency 
vector of Uq. Since b is constant and — > oo, we have limg_>.oo voi(rr/+2b) ^(^9) ~ P- 
Cq C n„(vol(mg + 2b) be the set of all possible color frequency vectors of all extensions 
of the colorings of Cr((mq),Cg) to Cr{{nig + 2b)). Clearly #Cg < (voi(m,+2b)+n-i^ ^ 

0(vol(mg + 2b)"^^), q ^ oo. For each c e n„(vol(m)), let Cr((m),c) be the set of those 
colorings in C'r((m)) that have color frequency vector c. From the above it follows that 

#Cr((m,),c,) < J2 #Cv{{mq + 2b), c) < #Cr((m, + 2b),£q)#C, 
cec, 

= #CT{{niq + 2b),eq)0(vol(mg + 2b)"-i) 

for some Cq £ Cq. Taking logarithms, dividing by vol(mg + 2b), and noting that 
Umg^.oo ™vo'i(m^^)'^^ = 1, we deduce that 

r Cg log#Cr((m, + 2b),^g) ^ ^ 

1™ — 77 — — ^TTT = P and lim -j^ — - — — — ^ — — > hrip) - e. 

q^oo vol(mg + 2b) q^oo vol(mg + 2b) 

Thus for p, q e Hr, e > we have sequences niq := (mi^g, . . . , md,q), := (ni,g, . . . , nd,g) € 
N'', 9 e N, with nig, ^ oo such that the following two conditions hold: 

Cr((m,),c,),Cr((n,),dg) 7^ 0,g e N, lim — rC,, = p, lim — J— = q, 

m,-)-oo vol(mq) ng^(Xi vol(nq) 

log#gr((m,),Cg) ^ , . . log#Cr(K),d,) 
hm I — — > hr{p) - e, lim -^^-^^ — ^ > /ir(q) - e. 

q^oc VOl(mg) q^oo VOl(ng) 

For i, j e N we show that 

hri^P + T^q) > ^hr{p) + T^/ir(q) - e- (4.13) 
t+j i+j i+j i+j 

Observe first that for any m, n e N** and c e n„(vol(m)) one has the inequality: 

#Cr((n.m),vol(n)c) > (#Cr((m), c))™'("). (4.14) 

Indeed, view as above, the box (n-m) as a disjoint union of vol(n) boxes (m). Color each box 
(m) in some color in the set Cr((m), c). Such a coloring is a member of Cr((n-m), vol(n)c). 
Hence (4.14) holds. 

Let ki,q := {mi^qni^q, . . . ,md^i,qnd~i,q),'kq := (ki^q, {i + j)'m.d,qnd,q). View (k,) com- 
posed of {i + j) boxes {niq ■ Uq). The above arguments show that 

#Cr((kg),vol(ki,g)(zc,+.7d,)> 

#(7r(((ki^q,imd,gnd,q), vol(ki,q)iCq)#C'r(((ki_,, jmd,gnd,,),vol(ki,,)jdg) > 
(#Cr((mg),Cg))^™'("')(#Cr((n,),d,))^™'(™'). 

Since Cr((kg), vol(ki^g)(iCg + jdq) D Cr((kq), vol(ki^g)(iCq + jdq) by considering the first 
coloring sequence in this inclusion, and using the maximal characterization of ^r(j^P + 

^q), we deduce (4.13). Since s was an arbitrary positive number we deduce (4.13) with 
e = 0. Since hr is upper semi-continuous we deduce the inequality hr{ap + (1 — a)q) > 
dhrip) + (1 — a)^r(q) for any a e [0, 1]. 

We now prove the claims (c-d). Assume first that Let u e diffPr- Then nr(u) = 
{VPr(u)} = dPr{u) and our assertions trivially hold. Assume next that the assumptions of 
part (e) of Theorem 4.3 hold. Recall that S{u) C nr(u) and conv S{u) = 5Pr(u) 2 nr(u). 
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Let Pi G ^(u), i = 1,. . So -Pr(u) = p^u + hr{pi),i = 1, • • ■ ,j- Since Ilr is convex, we 
obtain that for any a = {ai, . . . ,aj) & Hj p := ^iPi ^ ^r. As hr concave we deduce 

3 

Priu) =^aipJu + hr{Pi) < p^u + /ir(p). 

i=l 

The maximal characterization (4.6) yields that -Pr(u) = p^u + /ir(p)- So p e nr(u) and 
^r(p) aihriPi). This proves (c-d). 

We now prove (e). Recall that p G 9Pr(R") if and only if p G dPr{u) for some 
u e M". Use part (c) of this Theorem and part (c) of Theorem 4.3 to deduce the equality 
^r(p) = — -Pr(p)- If nr consists of one point then Hr = 9Pr(M") and (c) trivially holds. 
Assume that Hr consists of more than one point. Since 9Pr(K") 3 ri (domPp) ^ 0, use 
the second part of (a) of this Theorem to deduce /ir(p) = — -Pr(p) ea,ch p e ri (Hr). 
Suppose that q e nr\ri (nr),p € ri (Hr). Let 

/(a) := -/ir(aq + (1 - a)p), g{a) := Pr{aci + (1 - a)p), for a G [0, 1]. 

Since aq+(l— a)p G ri (Hr) for a G [0, 1) it follows that /(a) = g{a) for a G [0, 1). Since Pp is 
a proper closed function, it is lower semi-continuous. Hence nr(q) = ^(l) < liminfa/ii g{a). 
Since <? is a convex function on [0, 1] it follows that liminfayii (7(a) = limayi g{a) < g{l)- 
Hence g{l) = lima^ii ^(a). Recall that hr is a concave upper semi-continuous on Hp. 
Hence —hr is a convex lower semi-continuous function on Hr- Hence — /ir(q) = /(I) = 
lim„^i /(a). Therefore /(I) = g{l), i.e. /ir(q) = Pi^(q). ■ 

Wc now list several facts which arc consequences of Theorem 4.3. Given u G M", then by 
(d) each p G nr(u), namely each p achieving the maximum in (4.6), is a possible density 
of the n colors in an allowable configuration from Cr(Z'') with the potential u. That is, 
the relative frequency of color i is equal to pi. For each u where Pr is differentiable, there 
exists a unique density of the n colors. Assume that Pr is not differentiable at u. Then dPp 
consists of more than one point. Let S{u) be defined as in (e). Since 9Pr(u) = conv 5'(u), 
S{u) consists of more than one point. Hence by (e) Hr(u) consists of more than one point, 
that is to say, there is more than one density for u. In this case u is called a point of phase 
transition, sometimes called a phase transition point of the first order. 

Proposition 4.9 Let e :={!,..., G R". Then for allteR 

Pr(u) = t + Pr{u-te). 
Proof Recall the definition of Zr(m, u) given in (2.1). Clearly c((/))^e = vol(m). Hence 
Zr(m,u)= e'^^'^)^"^ ^ gc(^)^te+c(^)^(„-te) 

<peCr{(m}) 0eCr((m)) 

= e*™'('")Zr(m,u-te), (4.15) 

which implies the proposition. ■ 

Thus to study Pp, we may restrict attention to those potentials u = (mi, . . . , m„)^ that 
satisfy u„ = 0. (I.e. we reduce the number of variables in the function Pr to n — 1.) We 
show that the same holds for dPr and VPr- For u G M", we use the notation 

u := (ui, . . . ,u„_i)'^ 

for the projection of u on the first n — 1 coordinates, and extend it naturally to sets U := 
{u : u G U} for U C M". In the other direction, for x G M"~^, we use the notation 

t(x) := (xi, . . . , Xn-i, 1- Xi Xn-l)^ 
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for the unique lifting of x to the hyperplane S„ := {x e K" : x^e = 1}, and 

to(x) := (a;i,...,a;„_i,0). 
A straightforward computation shows that 

q^(z-z„e) = t(q)^z-2„ VqeM"-\zeM". (4.16) 
We define the convex function -Pr(') on R""^ by 

Pr(x) := Pr(to(x)) x e M"-^ (4.17) 
By taking i = u„ in Proposition 4.9, we obtain 

Pr(u) = M„ + Pr(u- u„e). (4.18) 

We now obtain a straightforward generalization of the density theorem proved in [12] for 

monomcr-dimcr tilings. 

Theorem 4.10 Let Pr{-) be defined on W'^ by (4.17). Then 

dPr{u) = dPr{u - M„e)) Vu e R", (4.19) 

aPr(K") = aPr(M""^). (4.20) 

Furthermore, Pr is differentiable at u if and only if Pr is differentiable atu — w„e. If Pr 
has all n — 1 partial derivatives atu — Une, then Pr is differentiable at u and 

(dP dP ^ ^ c) P I 

- • • • ' a^lfr^" " ' - § - ""^^ j ■ (^-'^^ 

Proof Assume that p e 9Pr(u). By definition of aPr(u) we have 

JV(u + v) > p^v + Pr(u) Vv. (4.22) 
Choose V such that Vn = 0. Then, using (4.18) for Pr(u + v) and Pr(u) in (4.22), we obtain 

Un + Pr(u - UnG + v) > p''"v + U„ + Pr(u - UnG) Vv, 

which by definition means p € aPr(u — u„e). 

Conversely, assume that q e aPr(u — Wnc), which means that 

Pr(u - u„e + v) > q'^v + Pr(u - u„e) Vv. (4.23) 

Now for arbitrary z, choose v = z — z„e in (4.23), and use (4.16) once and (4.18) twice in 
the resulting inequality to obtain 

Pr(u + z) - M„ - z„ > i(q)^z - z„ + Pr(u) - u„ Vz, 

which means by definition that t(q) G 5Pr(u), and therefore q S 5Pr(u). We have 
proved (4.19). 

We now show (4.20). It follows easily from (4.19) that aPr(R") C aPr(]R""^). To show 
the reverse inclusion, let q e 5Pr(w), and let u = io(w) so that w = u — Un^. By (4.19) 
we have q e 9Pr(u), as required. 

Assume that Pr is differentiable at u. Then aPr(u) is a singleton, and therefore 9Pr(u) 
is a singleton. By (4.19) dPriyi— Unc)) is a singleton, and therefore Pr is differentiable at 
u — u„e. 
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Conversely assume that Pp is differentiable at u — u„e, and therefore 9Pr(u — Une) 
is a singleton {q}. We assert that 9Pr(u) = {t(q)}. Indeed, if x e dPr{u), then x S 
dPr{u) = dPriu — u„e) = {q}, where the first equality is by (4.19), and so x = q. Since 
X S dPr{u.) C n„ C S„ by (f) of Theorem 4.3, it follows that x = proving the assertion 
since 9Pr(u) 0. Therefore Pr is differentiable at u. 

The last statement of the theorem follows from the previous statement and the fact [26, 
Thm 25.2] that a convex function / : M*^ ^ R is differentiable at a point a if merely the k 
partial derivatives i = l,...,k exist at a. ■ 

We can reformulate Theorem 4.3 for Pr. 

Theorem 4.11 Let Pp be the conjugate convex function of Pt- Then 
(a) /ir(p) < -%(p) for all p G Ilr. 
(b) 

Pr(x) = max(p^x+ /ir(p)) for all x e M""^ (4.24) 

pSlIr 

For X e M"~^, let q(x) e nr((-o(x)) i.e., q(x) is any vector satisfying 

t(q(x)) e Hr and Pt{^) = t(q(x))^to(x) + /ir(t(q(x))) 

= q(x) x + /ir(i-(q(x))). 

rc; /..r(4q(x))) = -p^(q(x)). 

fdj q(x) e 9Pr(x). /n particular, if x G diffPr, then q(x) = VPr(x). Therefore 
aPr(difrPr) C n^. 

('ej Lei x e R"~^ \ diff Pr, and let S{x) consist of all the limits of sequences VPr(xi) 
such that Xj e diff Pr and Xj ^ x. Then S{x) C nr(to(x)). 

(f) conv Ilr = domPp. 

Proof From (4.6) we have Pr(u) > p^u + hr{p) for all u G M" and p e Ilr- Fix p e Ilr, 
and let x e M"-^ and u = to(x). Then Pr(x) = Pr(u) > p'^x + /ir(p), so -/ir(p) > 
px— Pr(x). Now take the supremum over x € M""^ to obtain — ft.r (p) > -Fr(p)> which is (a). 
Substitute u = to(x) in (4.6) to deduce (4.24). We now show (c). By (4.25) and the definition 
of P* it follows that -/ir(t(q(x))) = q(x)^x-Pr(x) < Pp(q(x)). Combining this with the 
opposite inequality (a), we deduce (c). To prove (d), let x, z e R""^. Since t(q(x)) e Ilr, 
(4.24) applied to x + z and (4.25) give Pr(x + z) > q(x)^(x + z) + /ir(t(q(x))) = q(x)"^z + 
Pr(x). This proves (d). To prove (e), let x e M"-^ \ diff Pr and let q e S'(x). Then there 
exists a sequence Xj e diff Pr such that Xj — )■ x and VPr(xj) — > q. By the "furthermore" 
part of Theorem 4.10 we have to(x,) G diffPr and io(x) G M" \ diff Pr. By (4.21) and 
the continuity of t we have VPr(io(xi)) = t(VPr(xi)) t(q)- This shows that ;,(q) G 
5(io(x)) C nr(io(x)), where the inclusion is by (e) of Theorem 4.3. It follows that q = 
t{q) G nr((.o(x)), as required. Now we prove (f). By (a) we have Pr(p) < —hr{p) < oo 
for all p G Ilr, so Ilr Q domPp, and taking the convex hull gives conv Ilr C domPp. It 
remains to show the reverse inclusion. We know by Lemma 1 that 9Pr(M"^^) is the set of all 
points where Pp is subdifferentiable, so in particular ri (domPp) C (9Pp(M"^^). By (4.20) 
and (f) of Theorem 4.3 we have dPr{M."-~'^) = dPriR") C conv Ilr = conv U^. Combining 
the above inclusions gives ri (doniPp) C conv Ilr. Applying the closure operator, we obtain 
domPp C cl (domPp) = cl (ri (domPp)) C conv Ilr as in the proof of (g) of Theorem 4.3.^ 
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Since a probability vector p G n„ is determined completely by its projection p on the 
first n — 1 components, we can view the function hr ■ Hr — ^ R+ as a function on Ilr- 
Formally, let 

/ir(q) := /ir(i(q)) for all q e 11^^. (4.26) 



5 Pr and density entropies for one dimensional SOFT 

In this section we apply the results of Section 4 to one dimensional SOFT. In this case F is 
given by a digraph F := Fi C (n) x (n). 

Theorem 5.1 Let F C (n) x (n) be a digraph on n vertices, with at least one strongly 
connected component. Then -Pr(u) — logp(£'(F, u)), where the nonnegative matrix D{T,\i) 
is given in Proposition 2.2. IfT is strongly connected, or more generally F has one connected 
component, then Pr is an analytic function on M", Hr is a closed convex set of probability 
vectors equal to domPp, hr is concave and continuous on ri Fir, and coincides there with 
—Pp. In particular, for any u e M" 

ftr(p(u)) = -p(u)^u + logp(D(F, u)), where p(u) := ^^^i^^^'^j^ - (5.1) 

p(£>(F,u)) 

Furthermore 

hr = max hr{p) = hri^^^^^f^^) = logp(F). (5.2) 

Assume that F has k > 1 connected components Ai, . . . , A^. Then Pr(u) = 
max(PAi (u), . . . , PAk (u)), where each Paj is an analytic function on R". 

Proof Proposition 2.2 yields the equality Pr(u) = logp(£)(F,u)). Assume first that F 

is strongly connected, which is equivalent to the assumption that the adjacency matrix 
D(T) = {dij)ij^(^n) is irreducible. From the definition of £'(F,u) = {dij •e^^'^i "')jjg^„) 
it follows that £)(F, u) is an irreducible matrix for each value u e M". Then p{D{r,\\)) 
is a simple root of the characteristic equation p(z,u) dct(z/ — _D(F,u)) = for each 
u e R". Since the coefficients of p{z,\i) are analytic in u, where u S C", the implicit 
function theorem implies that p{D{T,\i)) is analytic function in a neighborhood of M" of 
C". Since p{D{V, u)) is positive on R" it follows that log p(_D(r, u)), u € R" has an analytic 
extension to some neighborhood of R" in C". Hence Pp is analytic on R". In particular, Pr 
is differentiable on R". Theorem 4.3 and Lemma 1 yield that domPp = conv Hp 2 Hr 2 
9Pr(R") 3 ri (domPp). Since Hp is closed we obtain that Hp = cl (domPp) = conv Hp, 
hence Fir is convex. According to Theorem 4.3 cl (domPp) = domPp. 

Since Pr is differentiable. Theorem 4.3 yields that /ir(p) = --Pr(p) for P e dPr{W). 
As 9Pr(R") 2 ri (domPp) we deduce that hr = -Pr on ri (domPp). Since Pp is a convex 
continuous function on ri (domPp), it follows that hr is a concave continuous function on 
ri (domP*). 

As Pr(u) = \ogp{D{r,u)) it follows that VPr(u) = Hence (5.1) holds. 

Clearly, p{D{r, 0)) = p{T) and (5.2) follows. 

Assume now that F is not strongly connected digraph. Rename the vertices of F such 
that D{r) is its normal form [13, Xni.4]. That is, -D(F) is a block lower triangular form 
matrix, where each submatrix on a diagonal block is either a nonzero irreducible matrix 
or 1 X 1 zero matrix. Then each nonzero irreducible submatrix corresponds to a strongly 
irreducible component of F. Let Ai, . . . , A^ be the k > 1 irreducible components of F. Since 
D(F, u) is also in its normal form it follows that p{D{r, u)) = maxjg[i^/j] p{D{Ak, u)). Note 
that logp(£'(A,,u)) = Pa.(u) for i = 1, . . . , fc. 

Assume first that k = 1 and Ai ^ F. Rename the vertices of F such that (to) is the set 
of vertices of Ai, where 1 < m < n. Let u = (wi, . . . ,Um)^ ■ Then Pr(u) = Pai(u) and the 
theorem follows in this case. 
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Assume finally that k > 1. The above arguments show that each Pa, is an analytic 
function in u, which does not depend on a variable uj if j is not a vertex of A^. ■ 

Assume that F is strongly connected and we want to compute p(u) = '^^D(r"u'))"' • ^® 
give the following simple formula for p(u) which is known to the experts. 

Proposition 5.2 Let T C (n) x (n) be a strongly connected digraph on n vertices. Let 
D{r, u) be the nonnegative matrix given in Proposition 2.2. Let x(u) = (a;i(u), . . . , a;„(u))^, 
y(u) = (yi(u), . . . , y„(u))^ be positive eigenvectors of D{T,u), D(r,u)^ respectively, nor- 
malized by the condition y(u)^x(u) 1. Then 

VPr(u) = = (yi(u)xi(u), . . . , 2/„(u)x„(u)) for each u e M". (5.3) 

Proof Let D(u) := D(r, u),p(u) := p{D(T,u)). Since p(u) > is a simple root of 
det(z/ — D{u)) it follows that one can choose x(u),y(u) to be analytic on M" in u. (For 
example first choose x(u),y(u) G R" to be the unique left and right eigenvectors of D{u) 
of length 1. Then let y(u) = ^^^^ j^^-^^ y(u).) Let di be the partial derivative with respect 
to Ui. Then 

y(u)^x(u) = IVu e R" ^ a^y(u)"^x(u) + y(u)'^aix(u) =0, f or i = 1, . . . , n. 

Observe next that y(u)^£'(u)x(u) = p{u). Taking the partial derivative with respect to Ui 
and using the formula (2.6) for the entries of D{u) we obtain 

d,p{u) = a,y(u)^i?(u)x(u) +y(u)^i?(u)a,x(u) +y(u)^9,;Z?(u)x(u). 
p(u)((9jy(u)^x(u) +y(u)^ajx(u)) + p{u)y,{u)x,{u) = p{u)y,{u)x^{u). 

This proves (5.3). ■ 

We now apply the above results to the following simple digraph on two vertices: 

Identify the red color with the state 1 and the blue color with the state 2, which is usually 
identified with the state 0. Then Cr(Z) consists of all coloring of the lattice Z in blue and 
red colors such that no two red colors are adjacent. This is the simplest hard core model in 

statistical mechanics. The adjacency matrix D{T) is the following 2x2 matrix ^ 



1 1 

Let u = {s,ty. Then VPr(u) = (pi(u),p2(u)) G Ha it follows that p2{u) = 1 -pi(u). 

It is enough to consider u — (s, 0) and pi{s) — '^^'"^lis'""' ■ So p :— pi{s) is the density of 
1 in all the configurations of infinite strings of 0,1, where no two 1 are adjacent. Clearly 
e* 



D{T,u) = ( j ^ . Hence 



, , 1 + VI + 4e^ , , 2e' 
PW = > Pi{s) 



1(1-^). (0,1, 



+ Ve-^ + 4)Ve-'* + 4 2^ y/l + Ae^^ 2' 



Note that pi{s) is increasing on M, and pi(— oo) — O,pi(oo) = i. Hence Hp — 

^ ri Hr. As Pr(0) 

(l + \/5)\/5 



conv ({(0,l)^,i(l,l)^}) and 9Pr(M^) ^ ri Hp. As Pr(0) = hv = log it follows 

that the value p* := pi(0) % ^ = .2763932024 is the density p* of I's for which 



hr^hriip*,l-p*))- 
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To find the formula for hr{p) = hr{{p, 1 — p)) first note that if p = pi{s) then 



Then 

/irW = logi^-plog^^i;^, pe(0,i). 

Our computations of Pr, for d>2, are based on upper and lower bounds, for example as 
given in Corollary 3.2. Wc claim that the function '°s^^^™^") (.^n be viewed as the pressure 
function of certain corresponding one dimensional subshift of finite type given. 

Consider for the simplicity of the exposition two dimensional SOFT given by F = 
(ri,r2), where Fi is a symmetric digraph. Let A be the transfer digraph induced by 
r2 between the allowable Fi coloring of the circle T(m). Then V := Cri,per('Ti) are the set 
of vertices of A. For any a, (3 G Cri,pei{'m) the directed edge {a, /3) is in A if and only if the 
configuration [(a, /3)] is an allowable configuration on Cr(('Ti, 2)). Note that the adjacency 
matrix D{A) = {dafi)a,0eCri,peT(m) is N x N matrix, where N := #Cri,per('Ti). Then the 
one dimensional SOFT is Cr(T(m) x Z): all F allowable coloring of the infinite torus in 
the direction 82 with the basis T(m). The pressure corresponding to this one dimensional 
SOFT is denoted by Pa(u). It is given by the following formula: Let 

^(A, U) = (rfaMu))a,/3eCp,.pe.(m), ^a^u) = d^0e'^(<^^+<^^^" ^ . (5.4) 

Then 

m 

The reason we divide logp(£>(A, u)) by m, is to have the normalization 

Pa(u + te) = Pa(u) + t for any t eR. 

It is straightforward to show that ^m"' "'' ^ ^a(u). Assume that A has one irreducible 
component. Then the arguments of the proof of Proposition 5.2 yield that -Pa(u) is analytic 
on R". Furthermore 

VPa(u) = (y(u)T(ai^(A,u))x(u),...,y(u)T(9„^(A,u))x(u)), (5.6) 

for any u G R". Here x(u) and y(u) are the nonnegative eigenvectors of £'(A,u) and 
D{A,u)^ , respectively, normalized by the condition y(u)^x(u) = 1. Then VPa(u) G Tin 
corresponds to the limiting densities of the n kind of particles in this one dimensional SOFT. 

In the numerical computations, as in the next section, we use one dimensional subshifts 
to estimate the pressure Pr from above or below as described for example in Corollary 
3.2. To estimate the partial derivatives of Pr one can find the partial derivatives of the 
pressure corresponding to the one dimensional subshift approximation using Proposition 
5.2. Since Pr(u) is convex in each variable we can estimate each partial derivative from 
above and below by finite differences. However, these estimates arc not as good as taking 
the derivatives of the one dimensional subshift approximation to Pr(u). 

6 The monomer-dimer model in 

A dimer is a union of two adjacent sites in the grid Z'^, and a monomer is a single site. By 
a tiling of a set 5 C Z'^ we mean a partition of S into monomers and dimers. By a cover 
of S we mean a tiling of a superset of 5* with each monomer contained in S and each dimer 
meeting S; in other words, dimers are allowed to protrude halfway out of S. Usually our 
set S will be a box or the entire Z**; in the case of a torus we only speak of tilings. As 
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mentioned in [12], the set of monomer-dimer tilings of 11^ can be encoded as an NNSOFT 
Cr(Z'') as follows. We color Z'' with the 2d+ 1 colors 1, . . . , 2(i+ 1: a dimer in the direction 
of efc occupying the adjacent sites i, i + is encoded by the color A; at i and the color fc + 
at i + e^; a monomer at i is encoded by the color 2(i + 1 at i. This imposes restrictions on 
the coloring, which are expressed by the d-digraph F = (Fi , . . . , Fd) on the set of vertices 
(2d+ 1), where 

(p, g) e Ffc ^ (p = fc, g = fc + d) or (p 7^ fc, g 7^ fc + d). (6.1) 

It is easy to check that this gives a bijection between the monomer-dimer tilings of and 
Cr(Z'^). Let Pr(u),u e M^'^+^ be the pressure function for the monomer-dimer model in 
. Since each dimer in the direction corresponds to the colors k and fc + it follows 
that fr(u) is effectively a function of d + 1 variables. To show that we define the following 
linear transformations 

Definition 6.1 (a) Let T,Ti : M^d+i the linear transformations 

T{wi,. . . ,Wd,Wd+i) = ( — ,...,—,—,...,— 
Ti{wi, . . . ,Wd,Wd+i) = {wi, . . . ,Wd,Wi, . . . ,Wd,Wd+i)- 

(h) Let Q : M^''+^ M'^+^ be the linear transformation given by Q{ui, . . . ,U2d+i) = 

(Ui + Ud+l, . . . ,Ud + U2d, U2d+l)- 

(c) Let Qd '■ K'' — > R &e the linear transformation (ui, . . . , Vdj^ i— > Wi + • • • + w^. 

Theorem 6.2 Let T = {Ti, ... ,T d)- coloring, with 2d + 1 colors given by (6.1) For 
u e let fr(u) denote the pressure function. Then 

(a) Pr{u)=PT{TQu). 

(b) aPr(R^'^+^) C r(R'^+i). 

(c) Hr = Tilid+i). Hence 9Pr(R'''+^) C T{I{d+i). 

(d) The function hr : Tlld+i M+ is a concave function. 

Proof (a) Since the colors i and i + d describe the two halves of a dimer in the direction 
ej for i = 1, . . . , rf, we have the identity Pr(u) = Pr{TQu). 

(b) Let p = {pi, . . . ,p2d+i) € 9(Pr(u)). In case u G diffPr, the Chain Rule applied to 

the identity in (a) yields the equalities pi = Pd+i, i = 1, . . . ,d. In case u G 9(Pr(u)) \ 
diff Pr, this follows from the fact that 3Pr(u) = conv S{u) as in the beginning of 
Section 4. 

(c) Since the color i appears with color i + d, it follows that pi = pi+d for i = 1, . . . ,d. 
Hence Ilr C Tlld+i. It is left to show that any 

p = {pi, . . . ,pd,pi, . ■ . ,Pd,P2d+i) e ^2d+i is in Ilr. Equivalently, the probability 
vector r := {2pi, . . . , 2pd,Pd+i) is the density vector of the dimer-monomer covering of 
U^. For d=l, this result is straightforward, e.g. [12]. So assume that d> 1. Suppose 
first that all the coordinates of r arc rational and positive: r = ( — , ^^i^), where 
m is a positive integer. Consider the sequence nig = {2qm, 2qm, . . . , 2qm) S N , g S N. 
Partition the cube (m,) to + 1 boxes with a basis (m^), (2gm, . . . ,2qm) G N"^"^: 
{{m'^,2qij)) , j = 1, . . . ,d + 1. Tile the boxes {{m'^,2qij)) with the dimcrs in the 
direction Bj for j = 1, . . . ,d, and the last box ((m^, 2qid+i)) with monomers. Then 
p = Tr e Ilr. Since Ilr is closed we deduce that Ilr 2 TUd+i- Hence Hr = TUd+i- 

In case u e diff Pr, then VPr(u) e Hr = TH^+i. In case u e a(Pr(u)) \ diff Pr, 
clearly 5(u) C Tlld+i. Hence 5Pr(u) = conv S'(u) C TUd+i. Therefore 5Pr(R2''+i) C 
T{Ud+i). 
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(d) According to part (b) of Example 4.7, the graph T, corresponding to the monomer- 
dimer model, is friendly, (as explained in [12, §4]). Part (b) of Theorem 4.8 yields that 
hr is concave Ilr = TUd- ■ 

Define 

iid(w) :=Pr(Ti(w)). 

In analogy with Proposition 4.1 and Proposition 4.9, i?d(w) : W^~^^ — )• M is convex Lipschitz 
function which satisfies the conditions 

|i?d(w + z)-i?4w)| < ||z|Uax w,zeIR'*+\ (6.2) 
Rd{w) = t + Rd{w - te) w eM.''+\,t eR. (6.3) 

We now derive the properties Rd{'w) that are analogous to the properties of i-r(u) 
discussed in Section 4. First we view Rd as the restriction of Pr to the {d + l)-dimensional 
subspace TR'^"'"^. Observe that QIl2d+i = QTIld+i = n^+i. Note that the vector r = 
{ri, . . . ,rd,rd+i)^ G H^+i can be defined intrinsically, where the dimer density in the 
direction e, for i = 1, . . . , d, and rd+i is the monomer density in the lattice Z'^. Let 

Hd{r) := hr(Tr), for any r e Ud+i- 

We view Hd{r) as the anisotropic dimer- monomer entropy of density r. 

It is straightforward to show that Rd satisfies an analogous theorem to Theorem 4.3. In 
particular, iid(w) = maxren<j+i(r''"w + Hd{r)). For w e we denote 

nd+i(w) := arg max (r^w + Hd{r)) = {r e n^+i : i?rf(w) = r^w + Hd{r)}. (6.4) 

Because of the equality (6.3), we can use the analogous results to Theorems 4.10 and 4.11. 
More precisely, for v G M"^, let Pd(v) be defined as in (4.17), i.e., 

PdM = Rdiioi^)) = Pt{{uo)), uJ = (v^,v^,0),v^ = (?;i,...,^;d)eR'^. (6.5) 

In other words, the two halves of a dimer in the direction of ej, are given the positive 
weight Xk — e"'' each, and a monomer is given the weight 1 = . Then Zpor(m,v) := 
Zr,per(m, Uo), is the grand partition monomer- dimer (counting) function in which we sum 
over all monomer-dimer tilings of the torus ^(m), and each tiling having exactly jii dimers 
in the direction for i = 1, . . . , d plus monomers contributes Jlf^j^ p^mvi ^ p^^ |--|^2], the 
function Z(m, v) := Z^ija., Uq) does not exactly count the weighted monomer-dimer covers 
of (m), because protruding dimers have only half of their weight counted. This can be easily 
taken care of as in [12], and the pressure Pd(v) is a convex function of v e M''. 

Lemma 2 Let v^ = (wi,...,Wd) G M'' and let cr : (d) — > (d) he a permutation. Then 
Pd{{vi, ■ ■ ■ iVd)^) = -Pd(('*^o-(i)! • • • )^^<7(d))^)) other words, Pd{^) is a symmetric function 
of v\, . . . ,Vd- Similarly for Z{in,v) □ 

Proof By applying the automorphism of N** given by 

(mi,...,TOd) 1-^ (TOa(i),---,TOa(d)) 

we obtain the equality 

Z{{mi, nid), (wi, . . . , Vd)^) = Z((to<^(i), . . . , m<^(d)), {v^(i),. . . , v„(^d)V), (6-6) 

and the result follows from (2.1). ■ 

Then each for each r e IId+i{Lo{v)) we have f G dPd{'v). We define := n^+i to be the 
projection of 11^+1 on the first d coordinates. Let 

hd{f) = Hd{r), rGHd+i. (6.7) 

We can repeat the proof Theorem 4.11 to obtain: 
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Theorem 6.3 Let be the conjugate convex function of the pressure function Pd- 

Then 

(a) ftd(q) < -Pdi^) for all q e A^. 

Pd(v) = max(q^v + /id(q)) for all v e W^. (6.8) 
For V e M'', we denote 

Ad(v) := arg niax(q'^v + hd{(\)), 

</iai zs to say q(v) S Ad(v) if and only if 

q(v)eAd and Pd(v) = q(v)^v + ;id(q(v)). (6.9) 

/i4q(v))=-P,*(q(v)). 

(d) Ad(v) C dPd{v). In particular, if v e diS Pd, then Ad(v) = {VPd(v)}. Therefore 
dPdidiS Pd) C Ad. 

(e) Let V e M*^ \ diff Pd, anrf let S{v) consist of all the limits of sequences VPd{vi) such 
that Wi e diff Pd and Vj v. Then 5(v) C Ad(v). 

(/J conv Ad = domPJ. 

Thus, the first order phase transition occurs at the points v where Pd is not differentiable. 

As in [15, 4, 17] we consider the total dimcr density q := qi + ■ ■ ■ + qd- This is equivalent 
to the equalities vi = ■ ■ ■ = Vd = v = log s, where s > is the weight of a half a dimer in 
any direction. We define preSd(w) := Pd((u, . . . ,v)^) = Pd(we) : K M. Then pres^ is a 
nondecreasing convex Lipschitz function satisfies |preSd(w) — \)veSd{v)\ < \u — v\. 

Proposition 6.4 For each e N Qd{^d) = [0, 1]. Let 

hd{p) '■= max hd{'C[), for each p [0,1]. (6.10) 

qeA<i,Qrfq=p 

Then 

preS(i(u) = laa^i pv + hd{p)- (6-11) 
pe[o,i] 

Furthermore, hd{p) is the p-dimer entropy as defined in [15] or [12]. 

Proof Let p G [0, 1] be the limit density of dimcrs, abbreviated here as p-dimer density, as 
(m) — ^ 00 as discussed in [12]. Wc recall the definition of the p-dimer density in terms of 
quantities defined in §4. (See in particular Definition 4.2.) 

For each m e N"^ and a nonnegative integer a G [0, vol(m)] define 

Cr((m), a) := Ue=(ci,...,c2d+i)en2a+i(voi(m))x2rf+i=oCr((m), c). (6.12) 

So Cr((m), a) is roughly equal to the set of all covering of the box (m) C Z'' with monomer- 
dimers, such that the number of monomers is a. (It may happen that some of the dimers 

protruding "out of the box (m), sec [12].) Then p €E [0, 1] is dimer density if there exists 
a sequence of boxes (nig) C N'' and a corresponding sequence of nonnegative integers Ug e 
[0, vol(mg)], such that 

m„ -> 00, Cr((m<,),a„) 9^ Vg G N, and lim = I - p. (6.13) 

9^-<x) vol(mg) 
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From the definition of the density set of the dimers it follows that p is a dinier density 
if and only p = Qd<\ for some q S A^. Since A,; = Hd+i, it follows that Qd{^d) = [0, !]• 
For each p e [0,1] let 

u(\ V log#Cr((m,),ag) 

hdVP) ■= sup limsup 77 1 ^ > 0, (6.14) 

m„a, 9-s-cx) VOl(mg) 

where the supremum is taken over all the sequences satisfying (6.13). Then hd{p) is the 
p-dimer entropy as defined in [12]. Let hpres^(p) be defined as in (6.10). We claim 

hd{p) = hd{p) for all p € [0, 1]. (6.15) 

Observe first that Cr((m),c) C Cr((m), C2d+i) for any c = (ci, . . . , C2rf+i) G n2rf+i(vol(m)). 
The definition of /id(q) and hd{p) implies straightforward the inequality /id(q) < /id(<9dq). 
Hence hd{p) < hd{p)- (6.12) yields the inequality 

nr^ n \ \ ^ /vol(m) + 2d\ 
#C7r((m),a) < ( ' jx 

max #Cr((m),c). 

C=(ci,...,C2(i+l)en2<i+l(vol(m)),C2(i+l=0 

Use the arguments of the proof of part (b) Theorem 4.3 to deduce the existence of q € Ag, 

such that Q^q = V and hd{p^ < ^d(q). Hence hd{p) < hd{p) and therefore hd{p) = hd{p)- 

To show (6.11), take v = i>e in (6.8). We have q^v = pv, where p = q^e. In (6.8) take 
the maximum in two stages. The first stage is for fixed p, and the second stage over all p.m 

The results of [17] yield that pres^(i;) is analytic. Since pres^(u) is also convex and not 
affinc it follows that pres^(f) can not be constant on any interval (a, &). Hence p{v) := 
pres^(t;) is increasing on M with p{—oo) = (no dimers) and p{oo) = 1 (only dimers). 
Therefore the analytic function p : M — ^ (0, 1) has an increasing analytic inverse v{p) : 
(0, 1) — )• K. Recall that pres2(p) is a convex function of p. Moreover 

dpres^ dv{p) dpres^ dv dv{p) dv(jp) 

__ = v{p)+p— —- = v{p)+p^-p^ = v{p). 

As v{p) is an increasing function oi p it follows that pres|J(p) is a strictly convex function 
on (0, 1). The corresponding dimer density entropy pres^(p) = —hd{p) is a strictly concave 
function. This is an improvement of the result of Hammersley [15] which showed that 
hd{p) is a concave function on (0,1). [11, Corollary 3.2] claims a stronger result, namely 
hd{p) + \{p\ogp+ (1 —p) log(l —p)) is a concave function on [0, 1]. (Observe that p\ogp + 
(1 — p) log(l — p) is a strict convex function on [0, 1].) 

Since pres^; is differentiable it follows that pres^(p) = pv{p) — preSj;(w(p)). Hence we 
obtain the well known formula, e.g. [4] 

hd{p{v)) = pres^(w) — p{v)v, where p{v) = pres^(w) for all w e M. (6.16) 

Note that hd{0) ■= limp^o hd{p) = and /id(l) := limp^i hd{p) is the d-dimensional 
dimer-entropy. 



7 Symmetric encoding of the monomer-dimer model 

The disadvantage of the encoding (6.1) is that the Ffe are not symmetric, so we cannot 
apply the results of Section 3 directly. However, as pointed out in [12], there is a hidden 
symmetry, which enables us to obtain results analogous to those of Section 3. We now 
adapt the arguments of [12, Section 6] to i-'d(v), the pressure corresponding to the weighted 
monomer-dimer coverings. 
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For d E N, K C {d) and m G N'^, we denote by (mx) the projection of (m) on 
the coordinates with indices in K. Let Cper,i<:(m) be the set of monomer- dimer covers of 
T(m/f) X (m^(i^\;f), and Zper,/f(m, v) the corresponding weighted sum. Thus Cper,(d)(m) = 
Cpcr(m) and Zp(,r,(d> (nij v) = Zpcr(ni,v). Note that by the isotropy of our T, ^^Cpoi-.K (m) 
is invariant under permutations of the components of m if if undergoes a corresponding 
change. Similarly for Zpor,i<-(m, v), if K and v undergo a corresponding change. 

In order to analyze Cpcr.{d}(Hi)j we focus on the dimers in the cover lying along the 
direction e^. More precisely, with m' = {mi, . . . ,md-i), we consider (m') x T(rnd) as 
consisting of md levels isomorphic to (m'). A subset S of the sites on level q is covered by 
dimers joining levels q — 1 and q (with level understood as level rrid); a subset T disjoint 
from S is covered by dimers joining levels q and g + 1 (with level rud + 1 understood as level 
1); and the remainder U of level q is covered by monomers and dimers lying entirely within 
level q. We are interested in counting the coverings of U subject to various restrictions. 
With that in mind, for m' G N''"^ we define an undirected graph G(m') whose vertices are 
the subsets of (m'), in which subsets S and T are adjacent if and only ii S f\T = $. When 
5 n T = we also define, using U = (m') \ (S U T), and v' = (ui, . . . , Vd-iY , 

asT^y') = sum of weighted monomer-dimer tilings of U 

^st(v') = sum of weighted monomer-dimer tilings of U viewed as a subset of 

r(m') 

Pst(v') = sum of weighted monomer-dimer covers of U, viewed as a subset of 
T(toi) X ((m2, . . . , rud-i)), each monomer within U, and each 

dimer meeting U but not SUT. 
cst(v') = sum of weighted monomer-dimer covers of U, each monomer within U, 

and each dimer meeting U but not SUT. 

In the tilings/covers counted by dsri'v'), &st(v'), Pst(v'), cst(v'), each monomer lies 
within U and each dimer meets U but not SUT. In asT(v'), each dimer occupies two sites 
of U that are adjacent in (m'). In &st(v'), each dimer occupies two sites of U that are 
adjacent in T(m'), so is allowed to "wrap around". In P5t(v'), the dimers in the direction 
of ei are allowed to "wrap around" and the other dimers are allowed to "protrude out" of 
((m2, . . . ,md-i))- In cst{'v'), the dimers may "protrude" out of (m'). The weight of each 
monomer-dimer cover is a product of the weights of dimers and "half" dimers appearing in 
the cover. If a dimer in the direction of is entirely within U, then its weight is e^'"''. If a 
dimer "protrudes out" in the direction of e^, then its weight is only 6"*=. Therefore 

asT(v') < 6st(v') < pst(v') < cst(v'). 

By definition, iiU = 0, then asri^') = ^st(v') = pst(v') = cst(v') = 1. Notice that when 
d = 2, there is no distinction between bsri^') and pst{^')- 

We define matrices A(m',v) = (asT(v))s,TC(m') > -S(m', v) = {bsT{^))s,TC{in'), P(m',v) = 
{psT{'v))s,TC{m') , C(m',v) = {csT{v))s,TC{m') with rows and columns indexed by subsets 
of (m') as follows: 
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-4(m',v)sT 



P{^m',^r)sT 



C(m',v)sT 



-B(m',v)sT 














if 5 n T = 
if 5 n T ^ 0. 



if 5 n T = 
if 5 n T ^ 



if S' n T ^ 
if 5 n T 7^ 



if 5 n T = 
if 5 n T ^ 



Thus ^(m', v). B(m',v), P(m'. v). C(m',v) are symmetric matrices — here is the "hidden 
symmetry" referred to above. Clearly 



(where the inequalities indicate componentwise comparisons). We use the notation a(m', v), 
/3(m',v), 7r(m',v), 7(111', v) for the spectral radii of these matrices, respectively, so that 



Note that by Kingman's theorem [21] all the spectral radii are log-convex in v. 

The four matrices have the same zero-nonzero pattern, namely the adjacency matrix of 
the graph G(m'). If the graph is connected, we say that the matrix is irreducible; if in 
addition the greatest common divisor of the lengths of all its cycles is 1, equivalently for 
sufficiently high powers of the matrix all entries are strictly positive, we say that the matrix 
is primitive. 

Proposition 7.1 Let2<deN and m = (m',md) G N''. Then 

(a) tr A(m', v)™'' is the sum of the weighted monomer-dimer tilings 0/ (m') x T{md) ; 

(b) trB(m',v)™<^ =Zpe,(m,v); 

(c) trP(m',v)™<^ = Zpe,,{i,d}(m,v); 

(d) trC(m',v)'"<^ = ^per,{d}(m,v) ; 

(e) fornid > 2, if column vector x{v) = {xs{v))sc{m') is given by xs{v) = bs${y')e'^^^'^ , 
then x(v)^B(m',v)™'«-2x(v) = Zper,{d-i)(m, v), 

if column vector y(v) = (ys(v))sc{m'> is given by ys{v) = cs${v')e*^'''^ , 
then y(v)''C(m',v)™<*-V(v) = Z{ni,v), 

and if column vector z(v) = (zs(v))5c(m'> is given by zs{y) = Ps0(v')e*'^''<*, then 
z(v)>(m',v)™--2z(v) = Zper,{i}(m,v); 

(f) the matrices A{in',v), B{m',v), P(m',v), C(m',v) are primitive. 

Proof We begin with proving (b), observing that (a), (c), (d) and (e) are similar. Assume 
first that nid = 1, and let (p € Cper(m). Since 4> can be extended periodically in the 
direction of with period 1, it can be viewed as an element of Cpor(m'). Therefore 
#C'per(m) = #Cpor(m') and moreover, Zper(m,v) — Zpor(m',v') (vd does not matter 
since no dimer lies in the direction of e^). We have tri?(m',v) = J2sc:(m') ^.s\s(v). Only 
the term 5 = contributes to the sum, and for this term we have U — (m') and = 



< A{in',v) < S(m',v) < P(m',v 



■) < C(m',v; 



a(m', V 



■) < /3(m',v) < 7r(m',v; 



■)<7(m',v). 
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Zpor(m',v') = Zpoi(m,v). Hence tri3(m',v) = Zpcr(m,v). Now assume that > 1, 
and consider a closed walk 81,82, ■ ■ ■ , Sm^jSi of length nid in G(m'). For each p' e 8q 
place a dimer in the direction of occupying the sites (p', q) and {p',q+ 1) (with rud + 1 
wrapping around to 1). We want to extend these dimcrs to a nionomcr-dimcr tiling of 
T(m') X Tinid) = T'(m), i.e., to a member of Cper(m), by monomers and by dimers not 
in the direction of e^, i.e., lying within the levels l,...,md- The weighted number of 
choices of such monomers and dimers to fill the remainder of level q is given by &sg_iS, (v'), 
and together with the weight of the dimers in the direction of intersecting level q it 
becomes 6sg_i5^(v). Therefore the weighted number of extensions to a member of Cper(m), 
i.e., the corresponding term of Zpor(m,v), is &SiS2(v)&S2S3(v) • • • &s„^_iS„^(v)&s„^Si(v). 
Conversely, each term of Zper(m, v) is obtained in this way. Hence Zper(m,v) is the sum 
of all the products of the above form, namely tr B(m', v)"*"*. 

To prove (f), we note that A(m',v) is irreducible, since whenever 5*0 T = 0, U can be 
tiled by monomers and therefore each subset of (m') is adjacent to in G(m'). Furthermore, 
A(m', v) is primitive since the graph has a cycle of length 1 from to 0. Since A(m', v) < 
S(m',v) < P(m',v) < C(m',v), it follows that B(m',v), P(m'.v) and C(m',v) are also 
primitive. ■ 

For the next lemma, we define Co (m) as the set of colorings of (m) corresponding to its 
monomer-dimer tilings (so no dimer protrudes out of (m)), and the corresponding weighted 
sum 

Zo(m,v)= e"^"^^^" = (v^,v^,0). 

</>eCo(m) 

Lemma 3 Let 2 < d & n and m! e M."^. Then 

log2o((m',md),v) . f , ^ /^ix 
lim =loga(m,v) (7.1) 

"id^oo rrid 

log^per,{d-i>((m',md),v) , , . 

lim =logp(m,v) (7.2) 

l0gZpe,,|i}((m',TOd),v) , , , . 

lim = log 7r(m , v) (7.3) 

m<i^-c!o rUd 

logZ((m',md),v) , / , N 
lim ^^^^ — ^=log7(m,v) (7.4) 

ma^oo md 

□ 

Proof From Part (a) of Proposition 7.1 Zo((m', rrid), v) < tr A(m', v)™"* , and therefore 

logZo((m^mrf),v) ^ logtr A(m^ v)'"-' 
lim sup ^-^^ — - < lim sup ^ — = log a(m , v). (7.5) 

rrid^oo md m^— >oo md 

The equality in (7.5) follows from a characterization of p{M) for a square matrix M > 0, 
namely p{M) = limsupfc_^3o(trM'^)fc (see for example Proposition 10.3 of [10]). Since 

— log Zo(m', md) is subadditive in rud, the first lim sup in (7.5) can be replaced by a lim. 
In order to prove the reverse inequality and thus (7.1), observe that each monomer-dimer 
tiling of (m') x T{md) extends to a monomer-dimer tiling in Co(m',md + 1) having the 
same weight (replace each dimer occupying (m', 1) and (m', md) by a monomer occupying 
(m', 1) and a dimer occupying (m', md) and (m', md -\- 1), and tile the rest with monomers). 
Hence Zo((m', -I- 1), v) > tr A(m', v)'"'* by Part (a) of Proposition 7.1. Therefore, since 

— logZo((m',md),v) is subadditive in md and thus the limits below exist, we obtain 

logZo((m^mrf),v) ^ ^.^ logZo((m^mrf + 1), v) 

md^oo md rrid^oo md 

logtrA(m^v)'"■^ 
> limsup = loga(m , v). 

TO(j— >(x) md 
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To prove (7.2), (7.3), (7.4), we use the fact mentioned in the proof of Proposition 2.2 that 
if M > and w is a column vector with positive entries, then p{M) = hmj;_>.oo(w M'^w)^. 
Applying this to M = B(m', v), P(m', v), C(m', v) and using Part (e) of Proposition 7.1 
with w = x(v),z(v),y(v) defined there proves (7.2), 7.3), (7.4). ■ 

Now we introduce the following notation. For m e N"^ and k G (d), 
jjj~fc := (toi, . . . , TTife-i, TTife+i, . . . , TTid) G N'^""'^. As special cases wc have the previous 
notation m' = m~'' and m~ = m~^. For v = {vi, . . . ,Vd)'^ G we use the notation 
v*^ := {vi, . . . ,Vk-i,Vk+i, ■ ■ ■ ,Vd,Vk)''' . Note that v** = v. Part (b) of Proposition 7.1 
implies 

Zper(m, v) = tr B(m', v)™'' = tr S(m^^ v'=)™'= . (7.6) 

Proposition 7.2 Let m e N'^,v G R'', and assume that is even. Then each k G 
{d — 1) satisfies 

log/3(m~'^,v) ^ log2 _^ log/3(m~'',v'=) 
vol(m~'') ~ ruk vol(m~'=) 

Proof We have 

/3(m'^'^,v)'"'' < trB(m'^'^,v)'"'' = tr B(m'~^'=, v*)"*' < 2™i('"~')/3(m'^'=, v*)'"^ 

The first inequality above follows since j5(m~'*,v) is one of the eigenvalues of S(m'"'^,v), 
which are all real, and rud is even; the next equality from (7.6); and the last inequality 
from the fact that i?(m~*^, v*^) has 2™^^™ ^ eigenvalues, all real, whose moduli are at most 
/3(m'"'^, v'^). Taking logarithms and dividing by vol(m), we deduce (7.7). ■ 



We define 



13 / \ V l0g^per,{l}((wi,m ),v) 

Pd_i(mi,v) := hm r , mi e N (7.8) 

vol(m-) 



m — >cx) 



Pd_i(0,v) :=log2. (7.9) 

Notice that for mi G N, Pd_i(mi,v) is the same as Pr(mi,u) defined in (3.1), where T is 
given by (6.1). For this reason the limit Pd_i(mi,v) exists. The following theorem is an 
analog of Theorem 3.1 and (3.11). 

Theorem 7.3 Let 2 < d G N, p,r G N, q G Z+, v G R'^. Then 

Pd-i{2r,w) > p . . > Prf-i(p + 2g,v)-Prf-i(2g,v) 

2r - "^^ ' - p ■ ^ ' ' 

Let m' = (mi, . . . , m^-i) G and assume that mi, ... , m^-i are even. Then 

, , log/3(m',v) 

Pdiy) < ^7, A 7.11 

vol(m') 

Proof Since #Co(m + 21) > #C(m) as explained in [12], it follows that 2'o(m + 21,v) > 
Z{m.,v). Hence, as in [12, formula (4.6) and (6.19)] and by Lamma 3, 

loga(m',v) log/3(m',v) log7(m',v) 
Pd(v) = hm = hm = hm ' • (7.12) 

m'->cx) VOl(m') m'->oo VOl(m') m'->oo VOl(m') 

First we prove (7.11). Let m' = (mi, . . . , Wd-i) G N''"-^, mi,...,md_i even, and let 
s = (si, . . . , Sd-i) G be arbitrary. Set 

mi = (si, . . . , Sd-i,mi), m.2 = {s2, ... ,Sd, nil, 1112), 

md_i = (sd,mi,...,md_i). 
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Note that (7.7) with fc = 1 states that 



log J(ni',v) ^ log 2 ^ log ,i(m .V"'-] 



vol(m') mi vol(m~) 

Using it — 1 times along with s = m'^, = nij, = 1113, etc., we obtain 

\ogP{s,{v2,V3,. . . ,Vd,viy) ^ log 2 _^ log/3(m]", (t;3, . . .,Vd,vi,V2V) ^ 



vol(s) Si vol(m^ ) 

log 2 log 2 log/3(mJ), (1-4, . . . , va, ('i, ('2, t's)""" 



Sl S2 vol(m2 ) 



< • • • 



d-1 



< 



log 2 log/3(m',v) 



s,- vol(m') 



Letting s — )• 00 and using (7.12) and Lemma 2 for the left-hand side, we deduce (7.11). 

We now demonstrate the lower bound in (7.10). Let m~ G N'^"-'^, s G N, g € Z+. Assume 
first that g G N. Since 7(m~,v-^) = p(C(m~,v^)) and C(m^,v^) is symmetric, it follows 
as in the arguments for (3.9) and by the analog of (7.6) for C(m~,v^) that 

^^"^ '''^ - trC(m-,vi)2« Zp,,,{i}(2g,m-),v) " ^'•'"'^ 

Taking logarithms, dividing by vol(m^), letting — > 00, and using (7.12), Lemma 2 and 
the definition of Pd-i{mi,v), we deduce the lower bound in (7.10) for the case q E N. If 
g = 0, we have to replace the denominators in (7.13) by tr / = 2™'^™ ^ , and the lower bound 
in (7.10) is verified by (7.9). 

We now prove the upper bound of (7.10). Let = (^2, • • • ,V(i-i,vi)''' . For each m' e 
N''"^ we have 

Jim'yf^ < tiCim'yf^ = Zpe,,{d}((m',2r),vi) = Zpe,,{i}((2r, m'), v), 

where the inequality above is true because the eigenvalues of the symmetric matrix C(m', v^) 
are real and 7(m', v^) is one of them, the first equality follows from Part (d) of Proposition 
7.1, and the last equality from (6.6). Therefore 

log7(m',v^) ^ logZper,{i}((2r, m'),v) 
vol(m') ~ 2rvol(m') ' 

and letting m' 00, we deduce the upper bound of (7.10) by (7.12), Lemma 2 and the 
definition of Pd- 1 (mi , v) . a 

In view of (7.12) we assume that '°vfi(™/'-)'^'* is a good approximation to Pd(v), and its 

partial derivative yoi(^jn')0{m' v) ^^^dv '^^ ^ good approximation to gj := ^^f^, the density 
of dimers in the direction of e,-. 



8 Numerical computations for the monomer-dimer model 
in 

In this section we explain in detail our computations for two dimensional pressure P2 (v) = 
P2{vi,V2) along the lines outlined in Sections 6-7. Our computations based on our ability 
to compute the spectral radius of the transfer matrix corresponding to the monomer dimer 
tiling of the torus (Z/m) x Z. This is a two dimensional integer lattice corresponding to a 
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circle of circumference m times the real line. In the notation of Section 7 this lattice is given 
by T(m) x Z. This transfer matrix is denoted by B{m,v). Let x = e^^,y = e^^. The the 
weight of the dimer in direction X, i.e. the horizontal dimer that lies entirely on the circle 
T(m), is .T^. The weight of the dimer in the direction Y, i.e. the vertical dimer that lies on 
two adjacent circles, is y"^. The matrix B{m, v) is of order 2™, corresponding to all subsets of 
(m). Denote by 2<™> the set of all subsets of (m). For S G 2<'"> denote by #5 the cardinality 
of the set S. Then B{m, v) = [y*'^+*'^ f{x, S, T)]s^Te20,.) . Here f{x, S*, T) = if 5 n T 7^ 0. 
For STlT = the function f{x, S, T) is a polynomial in x, which is the sum of the following 
monomials. Consider the set F := {m)\S U T viewed as a subset of the torus T{m). Let J" 
be a tiling of F with monomers and dimers. A dimer [i, i + 1], occupying spaces i,i + 1, can 
be in J^, if and only if i and i+1 are in F, where m and and m + 1 = 1 are adjacent. To each 
tiling corresponds a monomial x^', where I is the number of dimers in the tiling T of F. 
Then /(.x, S,T) is the sum of all monomials corresponding to all tilings of F. Note that if 
SnT = and SUT = (m) then f{x, S, T) = 1. Furthermore f{x, S, T) = f{x, T, S). Hence 
the transfer matrix B{m,v) is a nonnegative symmetric matrix. The quantity Pi(m, v), 
defined by (7.8), is given as the logarithm of the spectral radius of i3(mi, v). In numerical 
computations, we view ^^^^''"^ as an approximation to the pressure -P2(v). More precisely, 
one has the upper and lower bounds on the pressure which are given by (7.10). 

As in [12], the matrix B{m^v) has an automorphism group of order 2m, obtained by 
rotating the discrete torus T{m) and reflecting it. Thus, to compute the spectral radius of 
B{m,-v), it is enough to compute the spectral radius of the nonnegative symmetric matrix 
B{m,'v) whose order is slightly higher than . See for details [12, Section 7]. [12, Table 
1, page 517] gives the dimensions of B{m,, v) for m = 4, . . . , 17. We were able to carry out 
some computations on a desk top computer up to m = 17. 

We first apply our techniques to examine the Baxter computations in [4]. Baxter com- 
putes essentially the values of the pressure pres2(w) := P2{v, v) and the corresponding density 
of the dimers p(v) := ^pres2(t'). Recall that the corresponding density entropy h2{p{v)) is 
given by pres2(w) — vp{v) (6.16). Note the following correspondence between the variables 
in [4] and our variables given in Section 6: 

5 = 6", - =e-''+P^^«2('')^ p=-. 

s 2 

The case s = v = 00 corresponds to the dimer tilings of Z^. In this case p = p{oo) = 1 and 
/i2(l) has a known closed formula due to Fisher [8] and Kasteleyn [20] 
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r=0 

As in [4] we consider the following 18 values of s 

s"^ = 0.02, 0.05, 0.10, 0.20, 0.30, 0.40, 0.50, 0.60, 0.80, 
1.00, 1.50, 2.00, 2.50, 3.00, 3.50, 4.00, 4.50, 5.00. 

Wc computed the upper and the lower bounds for preSj (log s) for the above values of 
s, using inequalities (7.10) for v = (logs, logs) and to = 2, ... , 17. In these computations 
we observed that the sequence -Pi(^'''(^°s^'^°s^)) jg decreasing for r = 1, . . . , 8. So our upper 
bound was given by ^i(i6.('°g«.'"g«)) for all 18 values of s. The lower bound was given 
by ^i(i6.(i°g«.i°g^))-^i(i4.(i°g«.iog«)) ^-1 ^ 0.02,..., 0.3 and by A (17, (logs, logs)) - 
Pi (16, (logs, logs)) for other values of s. 

The values of Baxter for the pressure were all but two values between the upper and the 
lower bounds. In the two exceptional values s~^ = 1.5, 2.0 Baxter's result were off by 1 in 
the last 10th digit. As in Baxter computations, the difference between the upper and lower 
bounds grows bigger as the value of s increases. That is, it is harder to compute the precise 



30 



value of the pressure and its derivative in configurations where the density of dimers is high. 
This points to the phase transition in the case where 1? is tiled by dimers only [3, p'133]. 
The pressure value for = 0.02 computed by Baxter has 8 values. Our upper and lower 
bounds give 4 digits of precision of the pressure. For the value s = 1.0 our computations 
confirm the first 9 digits of 10 digit Baxter computation. (This value of the pressure is equal 
to the monomer-dimer entropy h2 discussed in [12].) For the values = 2.0, . . . ,5.0 our 
computations gives at least 12 digits of the pressure. 

We also computed the approximate value of the dimer density p(logs) = pres2(logs) 
using the following two methods. The first approximation was obtained by computing 
the exact derivative of -Pi(™'('°g''''°g-'^)) for rn = 2, . . . , 14. The second approximation was 
obtained by computing the raTio A(m,(iog(.+0,iog(s+t)))-Pi(r«,(iogs,iogs)) ^ ^ -^q_5 
m = 2, ... ,14. It turned out that the values of the numerical derivatives for m = 14 
agrees with most values of Baxter computations up to 5 digits, while the values of the exact 
derivatives agrees only up 2 digits with Baxter computations. Note that to compute the 
value of /i2(j3(log s)) we need the values of pres2(logs) andp(logs) (6.16). 

We next computed the approximate values of the pressure pres2((i'i, ■y2)) and its partial 
numerical derivatives for 18^ = 324 values. The 18 values of vi and V2 were chosen in 
the interval (—1.61,4.). (These values correspond to the 18 values of logs considered by 
Baxter.) For the lower bound and upper bounds we chose the values of 

M^M^^M)-m2M,.^ and MIl^ respectively. (8.1) 

Follows below the graph of -PiC^^'^Cji''"^)) ^nd the approximate values of h2{{pi,P2)), where 
Pi , p2 are the densities of the dimers in the direction xi , X2 respectively. The approximate 
values of /12 obtained by using the formula 

h2iiPUP2)) « ^^iHl^iiM -pi^i -P2V2, (8.2) 

Pi(14.(ri + /,.r,JJ- A(14.(ri.r.)) 
PI = Ylt ' ^^-^^ 

Pi{U,{vi,V2 + t))-Pi{U,{vi,V2)) 

= 14^ ' 

In our computation t = 10~*. For more detailed graphs with 42^ = 1764 points see 
http://www2.math.uic.edu/^friedlan/Pressurel7Jun09.pdf 

The graph of the pressure P2((a;i,a;2)) is convex and the graph of the density entropy 
h2{xi,X2) is concave. Both graphs look is symmetric with respect to the line Xi = X2. In 
reality this is not the case, since Pi (m, wi , W2) is the pressure of an infinite torus with a basis 
m. So in direction xi we have at most [^J dimers, while in the direction X2 we can have 
an infinite number of dimers. For m > 10 the difference -^i ("'■"i ."2)- Pi im,v2,vi) j^^^ than 

— m 

10~^, which explains the symmetry of our graphs. Note that in Figure 2 the densities pi,P2 
satisfy the condition pi,p2 G [0, l],pi +P2 € [0, 1]. The entropy h2 is in the interval [0, 0.67]. 

We also got similar graphs for the lower bound given in (8.1) and the corresponding 
analog of the approximation of h2{{pi,P2)) given by (8.2). These graphs were very similar 
to the graphs of ^^^^^'^^^''"^^^ and the approximation of h2{{pi,P2)) given by (8.2). 
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